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The cosmological gamma-ray burst (GRB) phenomenon is reviewed. The broad ob- 
servational facts and empirical phenomenological relations of the GRB prompt emission 
and afterglow are outlined. A well-tested, successful fireball shock model is introduced 
in a pedagogical manner. Several important uncertainties in the current understanding 
of the phenomenon are reviewed, and prospects of how future experiments and extensive 
observational and theoretical efforts may address these problems are discussed. 
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1. Introduction 

The study of astrophysical objects is overwlielmingly done by astronomers using 
the temporal and spectral information contained in the electromagnetic signals that 
these objects emit. Gamma-ray bursts (GRBs) are, by definition, electromagnetic 
signals in the gamma-ray band (in the spectral domain) with short durations (in 
the temporal domain). They are, however, unusual in having most of their electro- 
magnetic output in gamma-rays, typically at sub-MeV energies, and having most 
of it concentrated into a brief episode, typically lasting tens of seconds. 

The road leading to an understanding the nature of these objects has been 
bumpy, mainly due to the limited information contained in these abrupt gamma-ray 
episodes. For comparison, almost all the other astrophysical objects are observed 
in a broader spectral band during a much longer observation time. Historically, 
the lack (until recently) of observational breakthroughs on GRBs had left ample 
freedom for modelers to play around, so that the cumulative list of models champi- 
oned to interpret GRBs has been more numerous than for any other astrophysical 
phenomenon. As in other fields in science, experimental (or observational) progress 
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eventually constrains and eliminates most predictive models. At a certain stage of 
the development of a field, it is useful to ask oneself questions such as what is it that 
we know, what is it that we still do not know, and how could we find the answers 
for the unknowns. This is the main emphasis of this r eview . There exist alrea dy a 
number of comprehensive reviews on G RB obseryation J^ theories of GRI^^^ 

and on some general or selected topic J^ I ' l ^ l '^ l '' I I I Our purpose here 

is to provide a review which differs from others in several aspects. For example, 
we summarize the current observational progress in an organized, itemized form, 
making a distinction between solid facts and empirical laws. The introduction of 
the standard theoretical model is pedagogical, the content of which overlaps other 
reviews. However, we also dedicate a lot of space to discuss the uncertainties in the 
current models and highlight some active debates that are going on in the com- 
munity. We also attempt to reflect in an unbiased manner the work from various 
groups, rather than focusing on our own work, although some "selection effects" 
must still exist due to our incomplete survey of the literature. 

The structure of the review is as follows. We start with a brief summary of 
the GRB research history as well as the unique role of GRB study in general as- 
trophysics. We then compile the observational facts and empirical laws O, and 
introduce the standard theoretical framework, i.e. the fireball shock model (Q. In 
21 we outline some outstanding issues in the GRB field, which comprises most of 
the hot topics (and sometimes hot debates) in the current GRB scene. These issues 
include whether the fireball is dominated by the conventional (hadronic) kinetic 
energy content or by the less-studied magnetic energy; whether the GRB prompt 
emission occurs before the fireball deceleration (internal) or at the deceleration ra- 
dius (external); whether the GRB jets are uniform or with a quasi- universal angular 
structure; whether the GRB progenitors give rise to a one-step or to a (long de- 
lay) two-step collapse, etc. After outlining these issues, we attempt to foresee in 
13 developments in the upcoming new era for GRB study, led by future missions 
such as Swift and GLAST, complemented by other space- and ground-based detec- 
tors for both electromagnetic and non-electromagnetic signals. We review some of 
the theoretical predictions related to these observational capabilities and discuss 
how future observational campaigns could constrain and possibly settle some of the 
outstanding issues listed in ^ 



1.1. Major Landmarks in the Study of GRB 

As a framework for the following discussions, we list below some of the most impor- 
tant events in the development of the understanding of GRB. There may be many 
other important events which are not included in this list, and we apologize for 
any omission. Most of the papers selected below are marked by their high citation 
counts in the ADS archives, an indication that the works are significant and widely 
accepted. For the theoretical papers, we only include those that survived the later 
observational data and correctly contributed to the current knowledge of the GRB 
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phenomenon, or those which, ahhough not fully tested yet, represent "popular" 
views and draw broad attention from the GRB community. The papers extend up 
to October 2003, when this review is completed. They include both observational 
and theoretical developments, which have mutually influenced, and have benefited 
from, each other. 

A list of key advances on the observational side is: 



• GRBs were discovered s erendip itouslv in the late f960s, and the data were re- 
ported several years latej ^ l ' ^ 

• The Compton Gamma-ray Observatory {CGRO) was launched in 1991, with the 
Burst and Transient Experiment (BATSE) on board; BATSE provided evidence 
for an isotropic spatial distribu tion of GRBs, giving significant support to a 
cosmological origin interpretatior^l 

• Two categories of GRBs, "long" {t^>2 s) and "short" , {t-y<2 s) were identifiec^^ 

• Systematic analyses of GRB spectral data indicate that a so-called "GRB func- 
tion" (or "Band function") fits reasonably well most of the GRB spectra (for 
both classes of bursts 

• Another detector on board of CGRO, the Energetic Gamma Ray Experiment 
Telescope (EGRET), detected some bursts in the hard gamma-ray band; one 
burst, GRB 940217, was detected to have long-lived GeV emission extending for 
1.5 hripl; 

• A signature consistent with cosmological time dilation was detecteci^ 

• In 1997, the Italian-Dutch satellite Beppo-SAX pinpointed the first GRB low 
energy afterglow (for GRB 970228) in the X-ray band^, facilitating optical^ 
and radic^^ detections; 

• The first measurments of the GRB redshifts were obtained (for GRB 97050^^ 
and GRB 971214 ' , giving a solid proof that GRBs are at cosmological distance s: 

• A likely GRB-supernova association (GRB 980425 vs. SN 1998bw) was discoverecl^SEl 

• A bright and prompt optical flas h and a radio flare were discovered to accompany 
the energetic burst GRB 99012^SOEIl 

• An achromatic steepening break in the afterglow lightcurves was found in several 
bursts, hinting that at least some GRB fireballs are likely to be collimatedK- -; 

• X-ray spectral featur es wit h moderate significance were discovered in several 
GRB X-ray afterglowJSESl 

• Assuming the lightcurve steepening is due to the jet break, a collection of the 
burst jet data revealed that GRBs have a quasi-standard e nergy r eservoi m, 

• A class of so-called "X-ray flashes" (XRFs) were identified"-'^'""', which appear 
to be closely related to GRB; 

• Prompt optical afterglows from two more GRBs, GRB 021 004 and GRB 021211, 
were detected by ground-based robotic telescopeJ^^^^^^ 



^Notice that all the afterglow-related breakthroughs are for the so-called "long" GRBs (see il2lfor 
definition). No afterglow from a "short" GRB has so far been firmly identified. 
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• The p rom pt gamma-ray emission was reported to be strongly polarized in GRB 
02120^ (cf RefsHnilSI); 

• The first clear GRB-SN association was identified: an unambiguous supernova 
light-curve bu mp and spectrum in the GRB 030329 's optical afterglow lightcurve 
was discovered ^ ■ ' 1 ^ ^ 

• A distinct high energy spectral component which appears to evolve separately 
from the usual sub-MeV component was identified in GRB 94101/S3 

• Broadband observations of GRB 030329 suggest the there might exist a m ulti- 
component jet structure in some bursts, the total energy output being standar j ^^ l ^'^ l 

On the theoretical side, some of the major steps in the development towards the 
contemporary "standard" model are listed below. This list does not include many 
neutron-star-based Galactic GRB models (see Ref. 50 for a review), many of which 
were proposed before the proofs of a cosmological origin of GRBs became solid. 

• A generic argument about the "compactness problem" was raised^ in 1 975 ; 

• A self-similar solution of the relativistic blast waves was found in 197^. Al- 
though this study was not aimed at GRBs, it laid a mathematical foundation for 
relativistic GRB blast wave models, and has been extensively applied since the 
discovery of GRB afterglows; 

• Some generic features of a gamma-ray fireball were discussed i n 197 

• A cosmological model of GRBs was formally proposed in 1986 I I I involv- 
ing neutron star - neutron star (NS-NS) mergers as the power source. A pure 
photon-pair fireball was found to expand relativistically, leaving a quasi-thermal 
photospheric emission; 

• It was found that by including a small amount of baryon contamination, the fire- 
ba ll th ermal energy is essentially converted into kinetic energy of the expanding 
gaJS3, 

• It was po inted o ut that the kinetic energy is re-converted to heat and radiation 
via shockJ^SinS This suggestion laid the foundation for the current fireball-shock 
models. The specific scenario proposed in R.efs. 1581 l59l involves the shock forming 
when the fireball is decelerated by the ambient interstellar medium (ISM) , which 
is known as the "external shock scenario" . Synchrotron radiation was proposed 
as the mechanism of GRB emission, noting that this would soften and fade as a 
power law in time. 

• Around the same time, several major GRB progenitor models were further de- 
veloped, or proposed: the NS-NS merger model was further studied in some 
detah^^^ a "failed" supernova model was proposed^, which is now known 
as the "collapsar" model; and a millisecond magnet ar (n eutron stars with surface 
magnetic field of order 10^^ G) model was proposed^, 

• Another fireball shock scenario involving collisions among the individual s hells in 
the fireball outflow, known as the "internal shock scenario" was proposecl^^, 
as a way to resolve the very short time variability problem. 

• The hydrodynamical evolution of a fireball was extensively modeled both numer- 
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ically and analyticalljl™nHEnmmEaZ3 

The detailed reverse and forward shock synchrotr on and inverse Compton spec- 
tral components of external shocks were 

exploredZlESI, 

indicating the presence 

of a reverse shock optical synchrotron component (now identified as responsible 
for prompt optical flashes) and a GeV component due to inverse Compton. 
Based on the synch rotron radiation assumption, GRB radio and optical after- 
glows were suggested^^^^. A detailed self-consistent afterglow calculatioiP^stress- 
ing also the X-ray afterglow was performed shortly before the observational dis- 
covery of afterglows. The predictions include the power-law decay of the optical 
lightcurves as well as prompt optical flashes arising from the reverse external 
shock during the early stage of the afterglows. All these were verified by the later 

observations^^ 

Prompted by the discovery of the first GRB afterg the fireball model 

was studied and tested extensively, and found to be successful in interpreting the 
observation&22JML82^ 

It was suggesteci^ that GRBs are likely originated from "hypernovae" , and are 
associated with star forming regions; 

The simplest standard afterglow model was set uiJSSIHSl and was further devel- 
oped by incorporating additional effects^"^ 201 

Stimulated by the increasing evidence that GRBs are associated with supernovae, 
the GRB progenitor models involving collapses of massive stars were fur ther 
developed. The "coUapsar model" was extensively studied numericallj'^i'^^ESl 
and analyticallj^'^ A variant of the massive-star-collapse model, invoking a 
delayed GRB with respect to the SN explosion, was proposed^, known as the 
"supranova" model; 

The scenario involving GRB prompt gamma-rays from internal shocks and an 
afterglo w from the external shock as a unified sequence requiring a joint fit was 
analyzed '^ '^I'^^ t In parallel to this, extensive efforts were r nade to model the G RB 
prompt emission both within the internal shock modelJ^ '^l'*^*'^^^ ! '*'*'^'*' ! ^*'^"^ ! '*^^^'^ ! and 
the external shock modekl04 105 106 107. 

The c onsequen ces of GRBs originating from coUimated jets was propos ed and 
studiecP^^^^, leading to an extensive campaign of theoretical jet modelin j'*^^^^ ! '*^-*^'*^ ! '*^-*^"^ ! ^'*'^^ ! 
Prompted by the detection of the optical fiashes associated with GRB 99012^^, 
GRB 021004^9 and GRB 21211^0 411^ the reverse shock emission was stud- 
ied in greater detaifTTWHlTi6|117|l Isfl 19|1^()|1 ^111^2112^1. The recipes to use 
early afterglow data to systematically infer fireball parameters were recently 
proposea^^O -"^^^ ^^^; 

The possibility of producing Fe X-ray lines in supranova and in coUapsar stellar 
funnel models was investigatecPaiill^^Ml^l^^li^^l 

Broadband modeling of GRB afterglows was carried out, leading to constraints 
on some unknown model parameters , 

Prompted by the finding of a standard energy reservoir in long GRbJSSI, it was 
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proposed that long GRBs may have a quasi-universal structured jet configuratior^^^^^^ 
1.2. The multi- disciplinary nature of the GRB field 

The GRB field is almost unique in astrophysics in its multi-disciplinary nature. 
Involving stellar-scale events located at cosmological distances, GRBs straddle the 
traditional distance scales, and are a high energy phenomenon which emits a broad- 
band electro-magnetic spectrum extending over at least fifteen decades, as well some 
possible non-electromagnetic signals (such as cosmic rays, neutrinos and gravita- 
tional waves). This makes the GRB field an intersection of many branches in as- 
trophysics. 

GRBs were discovered at about the same time when quasars (or broadly speak- 
ing active galactic nuclei, AGNs) and pulsars, two other major developments in 
the 1960's, were discovered. The nature of those two classes of objects was agreed 
upon soon after their discoveries. AGNs are extra-galactic sources believed to be 
powered by gigantic black holes, while pulsars are compact neutron stars located 
in our Galaxy. Since the lack of GRB observational breakthroughs hampered the 
progress in this latter field, scientists from both the AGN and the pulsar commu- 
nity brought to bear their collective wisdom from their own fields to tackle the 
GRB problem. The GRB neutron star models were developed to fairly sophisti- 
cated levelJSni especially motivated by the reported detection of absorption and 
emission features in some GRB spectra-*^"^^ 1S9 turned out that the "classical" 
GRBs are of cosmological origin, so that the wisdom borrowed from the AGN com- 
munity (especially about blazars, the most energetic type of AGNs) finally bore the 
most rewarding fruit. Nonetheless, neutron star models eventually turned out to be 
partially c orrec t, i.e., one sub-class of such models that invokes ultra-strong mag- 
netic fieldJSni 

turned out to be useful in interpreting the class of soft, repeating 
GRBs (which appear to observationally distinct from the "classical" cosmologi- 
cal GRBs). These sources, known as "soft gamma-ray repeaters" (SGRs), together 
with another type of objects called "anomalous X-ray pulsars" (AXPs), are now 
believed to be "magnet ars" , neutron stars with super-strong (~ lO^'"^ G at surface) 
magnetic fieldJ^^^^^. Even within the current paradigm, pulsar models can still 
provide helpful insights for the understanding of cosmological GRBs. For exam- 
ple, the recent claim of strong polarization of gamma-ray prompt emission in GRB 
02120^ (cf. Refs. 031 144|l suggests a strongly magnetized central engine. A sim- 
ilar conclusion of a strong but less extreme magnetization was rea ched through a 
combi ned rever se-forward shock emission analysis for GRB 99012^^^ and GRB 

02i2iiinsinni. 

This hints that magnetic fields and/or Foynting flux may play an 
essential role in GRBs, and that some insights from pulsar wind nebula theories 
may turn out to be important to unveil the GRB mystery. We discuss this in more 
detail in 

Besides its close intimacy with the AGN and the pulsar fields, the GRB field is 
also broadly related to several other branches of astrophysics. First, in the stellar 
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context, since long GRBs are related to the deaths of massive stars, as supported 
by many observational fa cts, G RB study is closel y related to the fields of stel- 
lar structure and evolutio n I supernova q ' ■ ' ^ and supernova remnants (e.g. 
the external shock invoked in G RB models is the relativistic version of the super- 
nova remnant shodPSSHHini) . 

Progenitor studies have stimulated stellar popu- 
lation studieJl^. Central engine studies, on the other hand, have greatly promoted 
studies about mechanisms for extracting energy from accretion disks or spinning 
black hole&l^l ^^^l£li£li£f Second, within the galac tic context, GRB af- 
terglow lightcurve^^^^^^^ and spectral featureJ^^"^^^ probe the properties 
(e.g. density profile and chemical abundance) of the ambient interstellar medium 
(ISM) or the prestellar wind. Studying the properties of the GRB host galaxies 
as well as the GRB locations within the host galaxies brings valuable informa- 
tion about the global s tar forming history of the universe and the nature of GRB 
progenitorJ™SlIini. 

Third, GRBs are also important objects within the cos- 
mological context. Not only can they play a role similar to AGNs in probing the 
observed low-redshif t universe, but t hey may also form and can be detectable at 
much higher redshift j^^ ^4 | lt)5 | lt)t) | 

Thus, detecting high-z GRBs would allow us 
to see into deeper and earlier epochs of the universe, and to probe how the uni- 
verse ends its dark age and gets re-ionize Fourth, in high energy astrophysics, 
the origin of ultra high energy cosmic rays (UHEC Rs) has r emained a mystery for 
decades. Among many other models, a GRB origir^^^Sl [g a leading model for 
the so-called "bottom-up models" . In this scenario, cosmic ray protons are believed 
to be acceler ated in GRB sh ocks. The same processes can also generate high en- 
ergy neutrino j I I ^ 1 ^ ' ^ 1 go that GRBs are currently regarded by many as the 
top potential high energy neutrino sources for large area neutrino telescopes being 
built. Finally, sev eral GRB progenitor scenarios are believed to generate gravita- 
tional wave signalJP^^^, and GRBs are also one of the major targets for large 
gravitational wave detectors. On balance, it would appear that hardly any other 
field has such a broad range of interactions with other branches of astrophysics. 



2. Observational Progress: facts &: empirical relations 

GRB observational facts have been extensively reviewed, and here we summarize 
them in an organized, itemized manner. We will try to separate relatively solid 
observational facts from the empirical statistical laws which are derived from these 
facts. We divide our discussion into prompt emission and afterglow emission. 



2.1. Prompt emission 
2.1.1. Solid facts 

The main characteristics of GRB prompt emission have been collected in Ref. ^ 
(and references cited in that paper). Below we list the solid facts in itemized form. 
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Fig. 1. Typical lightcurve of a BATSE GRB, showing photon count as a function of time (in unit 
of second), from Ref. 1. 



including also some recent discoveries in the field''. 

• Temporal properties: 

(a) Durations (T, technically defined according to different criteria, e.g. Tgo or 
T50 being the time interval within which 90% or 50% of the burst fluence is 
detected): they span 5 orders of magnitude, i.e. from >,10~^ s to 10"^ s, typical 
values: '--^ 20 s for long bursts and ^ 0.2 s for short bursts; 

(b) Lightcurves: very irregular. Some bursts consist of very erratic, spiky compo- 
nents, while others are smooth with one or a few components. Some bursts 
contain distinct, well-separated emission episodes. Figure ^ is a typical GRB 
lightcurve. A compilation of different types of burst profiles can be found in 
RefCl 

(c) Widths of individual pulses (St) vary in a wide range. The shortest spikes 
have millisecond or even sub- millisecond widths, and St/T could reach as low 
as 10-3 - 10-^ 

(d) The vast majority of individual pulses are asymmetric, with leading edges 
steeper than the trailing edges, although only a small fraction can be visu- 



'^We discuss here only the prompt gamma-ray (and X-ray) emission. Prompt optical flashes have 
been detected in some bursts accompanied with the prompt gamma-ray emission. However, since 
they are generally interpreted as due to emission from the external reverse shock, we discuss the 
optical flashes in the early afterglow sub-section below. 
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Fig. 2. Examples for the GRB prompt emission spectrum (GRB 980329 and GRB 990510, from 
Ref.215). 



ally descerned. By in tegral means GRBs are asymmetric on all timescale^^^, 
and on averag e ! I Smooth single peak bursts typically have a fast-rise- 
exponential-decay (FRED)-type lightcurve. This comprises about 7% of all 
burstJIZ3. 

Spectral properties: 

(a) The continuum spectrum is non-thermal. Thermal (Planck-like) spectra are 
ruled out for the great majority of bursts. For most, the spectrum is well de- 
scribed by a smoothly-joining broken power law, known as a "Band-function'^^ 
(see Fig.El). Three independent spectral parameters are involved, i.e., a low 
energy photon spectral index (a), a high energy photon spectral index (/3), 
and the transition energy {Eq) or peak of the energy spectrum for (3 < —2 
(Ep). This spectral shape is valid both for the integrated emission over the 
whole burst duration, and for the emission during a certain temporal segment 
of the burst; 

(b) For a sample of 156 bright BATSE bursts with 5500 total spectreP^, it is 
found that a ^ —1 ± 1, (3 ^ ~'^-2^ the Ep distribution is lognormal, 
centered around ^ 250 keV with a full-width at half-maximum less than a 
decade in energy. The distributions of a and /3 are also generally suitable for 
describing fainter and softer bursts. The "narrow" Ep distribution among dif- 
ferent bursts, however, is likely to be influenced by selection effects. Va rious 
investigations indicate that the lack of high Ep bursts is likely intrinsicP"^. 
In the low energy regime, however, the narrowness of the distribution func- 
tion is mainly due to the "bright" flux-truncation in the sample. Lately, a 
group of X-ray transient events, which resemble normal GRB s in ma ny as- 
pects but with Ep around or below 40 keV has been 

identifieciSZEi. These 

bursts, named "X-ray flashes" (XRFs), extend the Ep distribution to the 
softer regime and are typically fainter. Globally, it seems that a "narrow" 
distribution of Ep may be obtained only when one is dealing with a sample 
of bursts within a narrow peak-flux rang&i^Q-lSi.. On the other hand, since 
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the bright BATSE sampl^^HI contains many spectra for each burst (at least 
eight), the results point to an intriguing conclusion, i.e., within a particular 
bur st, the JSn distribution is indeed very narrow. Also bursts tend to soft in 
see more discussions in 

(c) High energy spectral components: Dozens of B ATSE GRBs have also been 
detect ed at higher energies, e.g. by EGRE'M^ and Solar Maximum Mission 

Most of these detections are consistent with a Band spectrum 
extended to high energies without further breaks. Recently, however, a dis- 
tinct high e nerg y component was reported in the time-dependent spectra of 
GRB 94101'('^. This component sticks out at the high energy end of the 
conventional sub-MeV component, and extends to ~ 200 MeV (due to in- 
strumental upper cut off) with a photon number index of -1 that can not be 
attributed to an extrapolation of the MeV spectrum. More intriguingly, this 
component does not soften in time together with the low energy component, 
hinting an independent origin, i.e., perhaps a different emission site. At even 
higher energies (TeV), an excess of events coincident in time and space with 
GRB 970417 has been reported by the Milagrito groujU^. In the tempo- 
ral domain, GeV emission was detected in GRB 940217 lasting 1.5 houra^. 
However, this may be categorized as a high energy afterglow rather than 
prompt emission. 

(d) Spectral features: Absorption and emission features in GRB prompt emi ssion 
spectra were reported by the Sovi et sa tellites Venera 11 and Venera 1^^^, 
and by the Japanese mission The significance of these features was 
in the 2— 3(t range, with one claimed at 4cr. However, these were not confirmed 
by the BATSE tearriM The only spectral feature in prompt emission reported 
in recent years was a 3.8 keV absorption feature in GRB 990705 detected 
by BeppoSAX. 

• Polarization properties: 

Another important piece of information for electromagnetic radiation is its polar- 
ization, which is difficult to measure, especially in the gamma-ray band. Recently, 
a breakthrough in this direction was reported with the Ramaty High Energy 
Spectroscopic Imager (RHESSI). It was found that the prompt emission of the 
bright GRB 021206 was strongly polarized, and the claimed polarization degree 
IS 80%±20%^2I (cf. R.efs. Bnil^ . Unfortunately, due to its close angular distance 
from the Sun (which actually facilitated the detection by RHESSI) prevented the 
detection of the optical afterglow. A typical radio afterglow was however detected 
for this sourcel^. 

• Global properties: 

(a) Angular distribution: Isotropic for all bursts, or for cither long or short bursts, 
respectiveljECHi. 

There might be a small group of the so-called "long-lag" 

burstjni 

whose distribution is not isotropic but follows the super-galactic 
plane. Also, a sub-group intermediate between long and short bursts may 
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show departures from isotropj^^ 

(b) Intensity distribution: Generally there are two ways to quantify the intensity 
distribution. One way is through the peak flux distribution (i.e. the so-called 
log TV - log Pi^ plot>U. This results in a -3/2 power law slope at high fluxes, 
as expected in a homogeneous Euclidean model, and a shallower distribution 
at lower flux regimes that deviates from this simple homogeneous model. A 
second way is to find the average value of V/Vma,x: where V and Vmax are 
the volume of space enclosed by the distance of the source and the maximum 
volume of space enclosed by the distance at which the source could be still 
detected, respectively. It is found that < y/Knax > is less than 1/2 (the 
value expected for the homogeneous Euclidean model) , in the BATSE datM 
Such an intensity distribution, along with the isotropic angular distr ibution, 
are completely consistent with GRBs being cosmological event^^^ 

(c) Event rate: BATSE detected GRBs at a rate of about 1 per day. By correcting 
for sky coverage and other factors, the actual event rate is ~ 600 per year. Av- 
eraging over the whole Hubble volume (for Hq 70 km s~^ Mpc~^), this cor- 
responds to an average birth rate of ~ 7.5 Gpc^"'^ yr^^ or ~ 0.4 Myr~^ galaxy" 
where an average galaxy number density Ug ^ 0.02 Mpc"'^ is adopted. The 
local GRB birth rate is presumably smaller due to the drastic decrease in 
the st ar fo rming rate at low redshifts, and a widely quoted local (z = 0) 
rate i ~ 0.5 Gpc ^ yr ^, or ^ 0.025 Myr ^ galaxy ^. We now know 
that there are related cosmic events such as X-ray flashes (XRFs, see next) 
that may increase the total event rate by a factor of 2 or 3. This would also 
increase the GRB-XRF birth rate by a factor 2 or 3. We also know that 
GRBs are very likely collimated (see ^2.2|l . The beaming factor is uncertain 
dependent on the possible jet str ucture Assuming uniform jets, the true- 
to-observed beaming correction iJSHES ^ 500 for GRBs. For structured 
jets, this number is generally smaller. For a quasi-universal Gaus sian-t ype 
structured jet which fits reasonably well the current GRB-XRF dati^, the 
true-to-observed beaming correction factor is 14 for the combined GRB-XRF 
sample. The real GRB-XRF birth rate (which corresponds to the birth rate 
of their progenitor) should be increased by a same factor. Swift will help to 
determine both the observed GRB-XRF event rate and the geometric cor- 
rection factor more precisely. 

• Taxonomy: 

(a) The duration distribution is clearly bimodal (Fig. leading to a classifica- 
tion of GRBs into two types, i.e. "long" bursts with Tqq > 2 s, and "short" 
bursts with Tgo < 2 Such a classification is confirmed in the "hardness" 
domain, with long bursts typically being soft and short bursts typically be- 
ing harc"'^^. The long and short categories roughly consist of 75% and 25% of 
the total GRB population. It is commonly speculated that the two types of 
GRBs may have different origins. All our current knowledge of GRB coun- 
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Fig. 3. The bimodal distribution of durations of the BATSE GRBs, from Ref.189. 



terparts (afterglows and host galaxies, etc) are for long bursts. The nature 
of short bursts is still a mystery. There have b een extensiv e searches for a 
third type of GRBs with intermediate duration j -^^^^ l 194 | 1 95 |^ ^-j^e case is 
not conclusive; 

(b) Based on temporal characteristics (lightcurves), it is difficult to categorize 
GRBs into well-defined sub- types. Nonetheless, some phenomenological classes 
have been suggestea"'^. These classifications are however not fundamental. The 
differences among different types are likely caused by different behaviors of 
the central engine, different emission sites, or different viewing angles, not 
necessarily reflecting different types of progenitor; 

(c) Based on the spectral har dness, recently the so-called "X-ray rich GRBs" 
and "X-ray flashes" (XRFs) I I have been widely discussed as forming an 
apparently new type of transient event with respect to the conventional GRBs 
(Fig. 21). Whether they are the product of intrinsically different mechanisn is or 
are simply a natural extension of GRBs towards softer and fainter regimeJ^-^ 
is still unclear, but evidence is mounting that XRFs and GRBs are related 
events; 

(d) The "long-lag" bursts such as GRB 980425/SN 1998bw may belong to a sub- 
type of long bursts at closer distances (associated with the super-galactic 
plane) 



2.1.2. Empirical laws 

Besides the above relatively solid observational facts, some "secondary" empirical 
relations have been derived through various statistical analyses. These relations are 
potentially useful, posing important constraints on the theoretical models. Some 
are even helpful for deriving some important but unknown parameters. Below is 
a non-exhaustive list. The items 1-4 are for the temporal information alone, while 
the items 5-8 include the combined temporal and spectral information. 

• Power density spectra (PDSs) are powerful tools to quantify GRB temporal char- 
acteristics. It is found that the averaged PDS for many long bursts follows a 
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Fig. 4. The i?p-fluence distribution of the conventional long GRBs as well as the recently iden- 
tified X-ray flashes, from Ref.38. 



power law of index -5/3 oyer almost two decades in Fourier frequency, with a 
break around - 1 

at the higher end. This indicates that variabilities 
shorter than ^ 1 s are smeared out, for reasons as yet unclear; 
By tracking individua l pulses, the pulse width distribution for long bursts is 
found to be lognorma^^E^, with the peak variability timescale Stp]^ ~ 1 s. 
In some GRBs, there are so-called "quiescent times" where no gamma-rays are 
emitted are identified, and the durations of such quiescent episodes are positively 
correlated with the duration of the emission episode immediately following the 
quiescencJ^^. By subtracting th e quiesce nt times, the intervals between pulses 
also have a lognormal distributior^^EHOl 

The temporal behavior can be also quantified thr ough a "y ariabilitv" parameter, 
defined by some appropriate statistical definitionJ^EESI An intriguing empir- 
ical relation is that the more variable bursts (with larger V parameters) tend 
to have higher intrinsic isotropic luminosities L (derived from afterglow spec- 
troscopic measurements), i.e. L ~ V^'^. Although such a correlation is derived 
from a small parent sample, it exhibits a "Cepheid"-like correlation and makes 
"variability" a possible distance indicator of GRBs (Fig. bottom left) ; 
The cumulative lightcurves can reveal the global energy release rate of the central 
engine. For mos t bu rsts, the cumulative lightcurve has a constant slope as a 
function of timdl^EI There are two smaller groups (comprising ~ 4.8% and ~ 
2.8% of the whole population studied) whose cumulative lightcurves inc rease with 
time more rapidly or more slowly than the linear increase, respectivelj^^^ 
In many bursts (although not exclusively) there is a clear tre nd of spe ctral evo- 
lution. There are two types of evolution. One is "hard-to-soft'j^^^^SI in which 
hard emission leads the soft emission, and another is "tracking'^^^, for which the 
spectral hardness tracks the intensity. In either case, the Ep (derived by fixing 
both a and /3) decay is found to be exponential with photon fluence, and within 
a burst, the decaying constant is invariant from pulse to pulsJ^^I 
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• For asymmetric pulses, the pulse peak times niigratetowards later times at lower 
energies, and the pulse width becomes wi For 6 GRBs with spec- 
troscopic redshift measurements, it is found that the spectral lag is correlated 
with the isotropic luminosity as L53 « 1.3 x (r/0.01 s)~^ '^^, making spectral lag 
another possible "Cepheid-like" luminosity indicatoj '^1 (Fig. El top); 

• The low energy spectral index a is found to correlate withthe E'p of time-resolved 
GRB spectra, although the trend is still controversial^CI^El ______ 

• The Ep's are found to be positively correlated with the isotropic luminositjl^Sl^Sl 
With the spectroscopic redshifts of 12 BeppoSAX bursts, a firm correlation, i.e., 
Ep (X -^rad^^^ (where iJiad is the total isotropic energ y radiated in the X-ray 
and gamma-ray range, i.e. 1 keV - 10 MeV) is f oun a I which may be taken as 
yet another "Cepheid-like" luminosity indicatoJ^^ (Fig- HI bottom ri ght). The 
relation is found to extend from the hard GRB to the soft XRF regime02I^El 
However, due to the small sample effect, the issue is not conclusive, yet. Further- 
more, GRB 980425 had an Ep ^ 70 keV and Eiso 10**^ ergs, which apparently 
does not fit into the simple correlation; 

• Another redshift indicator is found in the gamma-ray spect ral d ata by taking 
into account of excesses from the exact power-law dependenc e ! '^ ^ ^\ 

• Luminosity function: Although the present afterglow spectroscopic redshift mea- 
surements provide too small a sample to perform a direct measurement, several 
attempts have been made to derive the luminosity function of GRBs. Although 
without full consensus, it is likely th at the luff ipositv is a power-law, with a 
possible steepening at high luminositie J''^'^''^l'^'^^^l'^'^^ l ^'''^ l '^^ 

• Using the V — L corr elation and the spectroscopic redshift dat JSSl, 

it is found 

that there might exist a cosmological evolution effect for GRB properties, such 
as luminosity (or even luminosity function), total energy, Ep^ duration, etc., al- 
though more spectroscopic redshift measurements are needed and selection effects 
need to be further clarifiecP^ESl. 



2.2. Afterglow emission 
2.2.1. Observational facts 

The afterglow observations have been extensively reviewed in Refs. El El El and 
references therein. Here we outline some of the main results in itemized forms. 
These include 

• Global properties: 

(a) Afterglows are (quite) broad-band, having been detected in the X-ray, the 
optical/infrared and the radio bands. In each band, the lightcurve generally 
displays a power-law decay behavior. X-ray afterglows are exclusively decay- 
ing when they are detected. Optical afterglows are generally decaying, with 
an initial early rising lightcurve having been caught in only a few bursts (e.g. 
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Fig. 5. Three "Cepheid"-like correlations found in GRB data, which could be adopted as rough 
distance indicators and redshift estimators. (1) Top: the spectral lag - luminosity correlation 
(from Ref.210); (2) Bottom left: the variability - luminosity correlation^"^ (from Ref.203); and (3) 
Bottom right: the Ep - luminosity correlation^ which extends to XRFs (from Ref.217). 

GRB 99012^. In contrast, an initial rising lightcurve in the radio band 
has been detected in many GRBs (followed by a canonical decay), an d the 
typical timescale for the radio afterglow to reach the peak is ~ 10 dayJ22l. 
Figure ini gives several examples of broadband afterglow lightcurves; 

(b) Very often there are various types of deviations from the simple power law 
decay. These include steepenings, bumps and wiggles (e.g. GRB 021004; GRB 
030329); 

(c) Not all bursts have afterglows detected in all of the three main bands. X-ray 
afterglows are the most commonly detected. About 60% of BEppoSAX bursts 
with X-ray afterglow detections are detected in the optical band. The other ^ 
40% of GRB afterglows are optically dark, i.e., there are no optical transients 
identified with sufficient exposure. These are dubbed as "dark burs ts" . In 
the HETE era, the fraction of dark bursts gets smaller, e.g. 109p2S. Radio 
afterglows are detected in about half of all GRB afterglows. 
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Fig. 6. Several well-monitored broad-band afterglow lightcurves. Upper left: GRB 970228 (Ref.79 
and references therein); Upper right: GRB 990510 (from Ref.33), notice the achromatic jet break; 
Bottom: GRB 030329 (from Ref.48), Notice irregularities in the lightcurves and the association of 
a supernova bump. 



(d) Essentially every GRB with an afterglow detection has an underlying host 
galaxy. The GRB host galaxy properties (magnitude, redshift distribution, 
morphologies, etc) are typical of normal, faint, star forming galaxies The 
GRB afterglow's positional offsets with respect to the host galaxy are con- 
sistent with GRBs being associated with the star forming regions in the 
galaxieJEOl 

(e) GRBs are at cosmological distances. Their redshifts are usually measured 
from the emission features of the host galaxies or the absorption features 
imposed on the afterglow continuum. As of October 2003, there are about 
33 redshift measurements, and detected redshifts range from 0.168 for GRB 
030329 (or 0.0085 for GRB 980425) to 4.5 for GRB 000131. See Refl^for 
a compilation of the afterglow data, including redshifts; 

(f) At least some GRBs are associated with supernova explosions. A famous 
example was GRB980425/SN1998bw association^^ 29 gjs^ igggbw was a pe- 
culiar, energetic Type Ib/c supernova. Using it as a template, other possible 
associations have been claimed by identifying a reddened bump in the op- 
tical afterglow lightcurves of GRB 980326^27^ GRB 97022&228.22l^ qrb 
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Fig. 7. The optical afterglow spectrum of GRB 030329 clearly reveals the spectral signature of 
Type-Ib/c supernovae (e.g. SN 1998bw) at day-10 since the burst trigger. This establishes a firm 
association between long GRBs and SNs (from Ref.45). 



000911^^ GRB 99120^^, GRB 9907li^, GRB 01112ll^, and GRB 
020403 I which may be attributed to supernova remnant contribution. 
Very recently, an unambiguous supernova signature has been detected in the 
z = 0.168 GRB 030329, firmly establishing the GRB-SN association in this 
objecl!S.46 (Fig. [71). 

• X-ray afterglows: 

(a) The continuu m sp ectrum is essentially a power law. By writing F^{t,v) cx 
t°'i'^ ^ one ha g I a ^ —0.9 and (3 ^ —1-4, both with a range of scatter. 
Late t ime lightcurve flattening was observed in a few bursts such as GRB 
00092^21 

(b) As of October 2003, X-ray emission line features have been claimed to be 
detected in the X-ray afterglows of 6 bur sts with mo derat e significance (< 
4.5cr, typically 3(t), i.e., GRB 970501*^, GRB 97082s3H3, GRB 99121^, 
GRB 000214233^ GRB 011211->^5., ^nd GRB 020813239; 

(c) Analyses of BATSE dat a rev eal soft tails in some GRBs which are likely to 
be early X-ray afterglowJ^^ 

• Optical afterglows: 

(a) Expressing the optical flux as i^opt(i, v) cx f^v^ , one haJ^^a ~ — 1 (at early 
times) and (3 ^ —0.7, both with a range of scatter; 

(b) In some bursts, a clear lightcurve steepening is seen after some time tbreak of 
order of days. The break is achromatic over different bands, and the temporal 
index after the break is typically —2 with some scatter. This break is typically 
attributed to the presence of a jet, and is termed the "jet break". At later 
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times, the decay rate gradually slows down until finally reaching a constant 
level due to the contribution of the host galaxy; __ 

(c) Other irregular temporal features are occasionally seen in various 

burstJ^m. 

These includ cP^^', for example, a substantial re-brightening in GRB 970508, 
an achromatic bump signature in GRB 000301C, wiggles in GRB 021004 and 
step-like features in GRB 030329; 

(d) Polarization: As of October 2003, 8 GRB optical afterglows have been de- 
tected to be polarizeG242 243 211^ j qj.^ 980425, GRB 990510, GRB 
990712, GRB 010222, GRB 020405, GRB 020813, GRB 021004 and GRB 
030329. The degree of polarization is small, typically several per cent. Large 
polarization angle changes (by nearly 90 degrees) were found in GRB 021004 
and GRB 030326^45 244. 

(e) Early optical flashes: Optical afterglow observations typically have started 
hours after the burst trigger due to technical reasons. However, early optical 
afterglows of two GRBs, i.e., GRB 990123^0 and GRB 021211^140^ were 
detected within 100 seconds after the GRB trigger. In both cases, the early 
afterglows indicate a steep decay with index ~ 2. For GRB 990123, a shajp 
rising lightcurve was also detected, and the V-magnitude is 9 at the peaM^. 
The early optical afterglow of a third burst, GRB 021004, was detected about 
3 mi nute s after the trigger, which displays a very shallow initial lightcurve 
decajEHl. 

• Radio afterglows: 

(a) In the radio band, the spectral index is generally positive for the observations 
which are typically in the 5 and 8.5 GHz bands. The lightcurves usually do 
not follow a simple power law declinJ^^ Some sources can be observed on 
timescales of years, and a late-time flattening (with respect to the standard 
fireball model) is often observea^47. 

(b) Prompt, short-lived radio flares have been detected in several burst J^SISHl 

(c) At early time s, rad io afterglows show strong fluctuations which are suppressed 
at later timea ^^. This can be interpreted as being due to interstellar scin- 
tillation effectJ^SI. The detection of such an effect in GRB 970508 clearly 
suggested that the source was expanding super-luminally, which gave a solid 
observational proof for the fireball shock modeP^. 

• Taxonomy: According to their afterglow data, long GRBs may be further classi- 
fied into several sub-types. Below are two such suggested classifications. They do 
not necessarily reflect intrinsically different groups (e.g. with different progeni- 
tors), but might be caused by environmental effects. So far there is no attempt to 
categorize GRB sub-types using combined prompt emission and afterglow data, 
but such efforts should be useful for identifying real sub-types. 

(a) Optically dark bursts: A fraction of GRBs with precise localizations do not 
have bright enough optical afterglows to be detectable. Possible reasons in- 
clude dust extinction, high redshift, or intrinsically faint nature. Evidence for 
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dust extinction has been collected for some burstJ^SHl Qn the other hand, 
observations for GRB 020124 and GRB 021211 indicate that at least someop- 
tically dark bursts are simply bursts with intrinsically faint afterg 
Indirect evidence for high rcdshift is also available for the recent dark burst 
GRB 031026255. 

(b) Fast-fading GRBs: Several bursts (e.g. GRB 980519, GRB 980326) show a 
steep afterglow decay (a ~ —2) in their early phase. They do not fit into the 
"standard energy reservoir" scenaricP^ and may constitute a peculiar class 
of GRBs256.. 



2.2.2. Empirical laws 

As in the prompt emission studies, some secondary empirical laws have been dis- 
covered in the afterglow data. Below is a non-exhaustive list. 

• Perhaps the most intriguing finding is that long GRBs seem to have a standard 
energy reservoir. This conclusion is based on the model that GRBs are coUimated 
and that lightcurve steepenings are due to a jet. From a simple toy model, one can 
derive the so-called jet opening angle 9j using the jet break time. For those bursts 
whose redshifts, and hence the total gamma-ray energies (assuming isotropic 
radiation) in the prompt phase, E^^iso, are measured, it is found that E^^isoO'j is 
essentially a constant, which means that the geometrically corrected gamma-ray 
energy for different bursts is essentially an invarianlP^ElSI (Fig. (S)); 

• The total GRB fireball energy should be Etot > E~^ + Ek, where E^ is the energy 
released as gamma-rays in the prompt phase, and Ek is the kinetic energy which 
is left over after the prompt phase and is dissipated in the afterglow phase. The 
lat ter conip onent (Ek) may be directly measured using broadband afterglow 
fit J ''^ I I or via analyses of the X-ray afterglow data for different GRBs at the 
same afterglow epoch257 258 259_ goth approaches lead to the same conclusion 
that Ek is also a standard value for different bursts. Combining this and the 
previous argument suggests that long GRBs have a standard energy reservoir; 

• An anti-correlation has been derived between the GRB spectral lag and the jet 
opening angle This relation, when coupled with the Cepheid-like correlation 
for the spectral laj^^, is another manifestation of the anti-correlation between 
the isotropic luminosity and the jet open ing angle, which is the direct consequence 
of the standard energy reservoir relation"^ 

• Other model parameters can also be inferred from broadband afterglow fits. Con- 
sensus is emerging about the values of some of these parameters. For example, 
the environment of most GRBs appears to be an interstellar medium (ISM) with 
an approximately constant density, although the afterglows of a small group of 
GRBs are consistent with a stratified wind-type medium (with a density profile 
n (X 7--2| ioi | io^ | xhe ISM density appears not to vary greatly among burst J I 
In the shock region, the portion of the energy which goes into the electrons (de- 
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Fig. 8. The geometry corrected gamma-ray energy (i.e. Ej 



?/2, where E-yjso is the 



total energy emitted in gamma-rays assuming isotropic radiation, and 9j is the jet opening angle 
inferred from afterglow lightcurves) is found to be a constant in many bursts, referring to a 
standard energy reservoir of long GRBs''^. Shown is the distribution of E-y with the latest data 
(from Ref.256). 

noted as £e) is typically 0.1 with some scatter, while the energy portion that 
goes into the magnetic fields (denoted as 63) is typically ~ 0.01 or less^"^^ ^23. 



3. Theoretical Progress: Standard fireball shock model 

In this section, we briefly introduce the key ingredients of a generic standard GRB 
fireball shock model in a pedagogical manner. This theoretical framework is the 
most widely used for interpreting the current GRB and especially afterglow obser- 
vations. Its "standard" nature is the product of its predictive power and success 
in passing various observational tests. This section focuses on those aspects of the 
model which are robust and have the least uncertainties. Discussion of some of the 
less certain aspects of the fireball model are deferred to 21 The contents of this 
section overlap to some significant degree several other reviews such as Refs. 0] El 

m 



3.1. Relativistic bulk motion 

A first ingredient of the standard model is that the emitting material responsible 
for GRBs and afterglows must be moving rclativistically. This is a consensus of 
all cosmological GRB models (even if these models are non-standard and differ 
considerably in technical details) and even for the old galactic halo models. The 
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arguments leading to this conclusion are straig 

htforwardSaSlESEi, 

as discussed 

by, e.g. Refs.01inilZlEl For a typical GRB gamma-ray fluence F ^ 10~^ erg cm~^ 
and distance D ^ 3 Gpc, the total isotropic gamma-ray energy released is typ- 
ically E = AnD^F ~ 10^^ ergs. Naively (without relativistic motion), the scale 
of the emission area is c5t = 3 x 10^ cm ((5t/10ms). Assuming that a fraction 
fp of photons is above the two-photon pair production (77 — > e+e^) threshold 
[€162(1 — COS012) > 2(meC^)^, where mg is the electron mass, ei and £2 are the 
energies of two photons, and 9i2 is the angle between the momenta of the two pho- 
tons], and using an approximate pair production cross section of the order of the 
Thomson cross section (Jt = 6.25 x 10~^^ cm^, the pair-production optical depth 
is huge, i.e. r^^ = fparFD^ / [cStfmeC^ - W^fp{F/lQ-^ erg cm-^){D/i Gpc)^ 
{5t/lQ ms)~^. Thus, the gamma-rays should have been attenuated in the source 
before traveling through the universe and reaching the earth. The only way to get 
rid of this apparent paradox is by invoking relativistic bulk motion, i.e., the GRB 
emitting region as a whole moves towards us observer with a high Lorentz factor. 

The relativistic motion eases this "compactness problem" in two ways. Suppose 
that the emitting region (e.g. an ejected shell of relativistic material) flies towards 
the observer with a bulk Lorentz factor T. First, the photon energy is blue-shifted 
by a factor of F, so that in the shell comoving frame, the bulk of gamma-rays as 
observed are actually X-rays. This greatly decreases the number of photons above 
the pair production threshold, i.e. fp drops by a factor of r^^"~^\ where a ~ 2 is 
the observed photon number spectral index, i.e., N{Eph)dEph oc E~j'^dEph. Here 
the factor (a—1) arises from the integration leading to the above-threshold photon 
numbers, and the factor 2 takes into account the con tributions of both the test 
photons and the target photons for pair productiorP^. A second effect introduced 
by considering relativistic motion is that the real physical scale of the emission 
region is T^cSt for an observed timescale of 5t (see tl3.2l below). So altogether the 
pair optical depth drops by a factor oF r2+2Q! ^ p6 p^j, typical bursts, F > 100 is 
required to have r-y^ < 1. This is only a rough estimate; more sophisticated analyses 
result in various lower limits of F in different bursts^^^ "^^^ . The typical 

F required to satisfy the observations is of order ^ 100. Hence, GRBs involve the 
fastest bulk motions known so far in the universe. 

There are other evidence of relativitic motion of the fireball. One is from the 
radio afterglow data that initially shows large interstellar scintillation but gets 
supressed later on. This presents a cle ar ev idence of superluminal expansion of 
the fireball caused by relativistic motior!^^. Another is from the analyses of early 
afterglow reverse shock data-'^-'^^ 12Q xhcsc analyses directly point to a large initial 
Lorentz factor of the fireball (sec i)5.1l for more discussions). 



'^Notice that the expressions of the optical depth in some of the previous reviewJ^'is incorreci '^^^ l 
given the above definition of a. 
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Fig. 9. The geometric configuration among the GRB central engine, relativistic emitting shells 
and the observer. 

3.2. Two reference frames, three timescales 

A typical GRB problem involves three major physical elements (Fig. a central 
engine, a relativistically moving shell (ejected by the central engine) which produces 
the GRB emission, and an observer. There are essentially only two inertial frames. 
One is the rest frame of both the central engine and the observer (aside from a 
cosmological redshift factor, which is small compared to special relativistic effects), 
and the other is the rest (comoving) frame of the relativistic shell or ejecta. The 
physical quantities (e.g. scale length and time) as viewed in the two inertial frames 
are different, and are related through special relativistic Lorentz transformations. 
For example, a length scale A' in the comoving frame is converted to A = A'/F 
along the shell's moving direction in the rest frame of the central engine. 

Similarly there are only two sets of clocks attached in both inertia frames, so 
that dt' — dt/T, where dt and dt' are the time intervals elapsed for the same pair of 
events in the central engine/observer frame and the comoving frame, respectively. 
However, in the GRB problem, there is a third relevant time scale involved. The 
complication comes from the propagation effect. In general, we can consider a shell 
emitter moving with a dimensionless speed /3, at an angle 9 with respect to the line 
of sight of the observer (Fig.Ejl. In the rest frame of the central engine/observer, the 
shell emits a first photon towards the observer at the time ti at the location A (the 
radius r), and emits a second photon towards the observer at time ^2 at the location 
B (the radius r + rfr), as recorded by clocks precisely adjusted in this inertia frame. 
The time interval for emitting these two successive photons is dt ~ t2 ^ ti dr/c. 
Let's assume that the distance between the observer and the location A is L. The 
first photon arrives at the observer at ti = ti +L/c, while the second photon arrives 
at the observer at ^2=^2 + {L/c — Pfidi) where /i = cos 6. These two times are also 
measured by clocks precisely adjusted in the same inertial frame. The time interval 
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for the observer to receive the two adjacent photon signals is simply 

dt^ {1- (3fi)di ~ di/2T^ = dr/{2T^c), (1) 

where in deriving the second part of the equation, we have assumed F = (1 — 
^2-j-i/2 ^ 2 and 9 1 (so that fi ~ 1). Notice that eq. 1^) is a pure propagation 
effect, which is also valid for the non-relativistic case. It is not noticeable in the 
Newtonian case, simply because dt ~ dt when /3 <C 1. In the relativistic regime, 
this effect becomes very important, and eq.(^ is one of the fundamental arguments 
used to solve the "compactness problem" as discussed above in H3.ll 

It is worth noticing that dt and dt are two different times in the same inertial 
frame, which describe two different pairs of events, dt describes the time for the 
shell to emit two photons, while dt describes the time for the observer to receive 
the same two photons. Since dt also describes the actual time of the shell behavior 
(e.g. when the shell moves to the location A or location B), the best way to refer 
to these two times is to call dt the "time in the rest frame of the central engine" 
or "the time in the fixed (or lab) frame" , while calling dt the "observer's time" . 
In the literature, the latter is sometimes called "the time in the observer's frame" , 
which is in principle right, but this is not the exact distinction between dt and dt. 
Sometimes in the literature dt is called "the time in the burster's frame" . This can 
be confusing to some readers, since the "burster" may be understood either as the 
central engine powering the burst (in which case the definition is correct) , or it may 
be understood as the flying shell which emits the burst (in which case the definition 
is incorrect). Also notice that in some articles, the dt discussed here is denoted as 
dt, while the dt discussed here is denoted as dT, dtobs or dt^. Since the observer's 
time is the most relevant one to describe the phenomenon, it is more convenient to 
define it as dt, the most straightforward notation. In the rest of this review, we will 
adhere to such a notation system. 

A variant of dt which is also widely discussed in GRB problems is the so-called 
"angular time" dtang- Instead of fixing 9 and varying r (or t), in some problems 
one needs to fix r (and hence t) but vary 9. This is relevant for the problem in 
which a shell-like emitter (with a fixed radius r) is illuminated instantaneously at 
the same engine-frame t, while the observer sees a long-duration emission caused 
by the time delay of radiation coming from higher latitudes. In a similar manner as 
for the derivation oidt, one dtang = {r/c)sa\9d9, or tang = (r/c)(l — cos 6*). For 
relativistic motion (F ~ 1/6* ^ 1), one has 9 ^ sai9 ^ 1/F, so that tang — (f/cF^), 
which is of the same order as t (eq.^)^*^*^ . Notice again that tang is also a pure 
propagation effect and is also applicable for the Newtonian case. 

The third time scale is the comoving time of the shell, measured by another set 
of clocks. It is related to the engine-frame time dt through 

dt' = di/T = dt/T{l - = Vdt ~ 2Tdt, (2) 

where V = [F(l — /3^)]^^ is the Doppler factor (since the observed radiation fre- 
quency is boosted by the same factor with respect to the frequency in the comoving 
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Fig. 10. Various characteristic radii in a generic relativistic fireball. 



frame, i.e., v — Vv'). The final approximation in eq.jSl is again for F 3> 1 and ^ ^ 1. 
3.3. Fireball evolution, characteristic radii and times 

The evolution of TOneric fireballs gxpanding i n^ an ambient medium has been 
extensively studiej^ ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^-'^ I ^^ ESl^ assuming a fireball composition of 
photons, electron/positron pairs and a small amount of baryons (but negligible 
magnetic fields). In a simplest toy model, we assume the following input parameters: 
(a) an average constant luminosity of the central engine, L; (b) the duration of the 
central engine energy injection, T, so that the total energy of the fireball is i? = LT 
and the initial width of the whole fireball shell in the fixed frame is Aq = cT; (c) the 
variability time scale, \ ms ^ T, which is due to the intermittent nature of the 
fireball central engine, as refiected by the spiky, irregular GRB lightcurves f H2.1.1|l . 
The fireball shell therefore actually consists of many mini-shells; (d) the average 
mass loading rate, M, so that the so-called dimensionless entropy \s rj = L/Mc^\ 
and (e) the particle (usually hydrogen) number density of the ambient medium, 
n (in the simple toy model here, n is taken as a constant, which is typical for an 
interstellar medium, ISM). The evolution is characterized by several characteristic 
radii (see Fig. 1101 for a cartoon picture), which we will discuss in turn below. 

• Initial state {r — tq): The base of the fireball fiow is connected to the GRB central 
engine, a black hole - torus system or a rapidly rotating magnetar. Ass umin g 
a IOM0 black hole, the initial length scale (taken 3 Schwarzschild radiij^ is 
To ~ QGM/c^ ~ 10^ cm. This radius is also the typical width of the mini-shells, 
i.e. (5o = cty = (3 X 10^ cm) {ty/1 ms). The initial state of the fireball is hot, 
with photons and pairs in equilibrium with a temperature To = (i/47rrQ(T)^/'* ~ 
(10^° K) Lll-^r'^y^ . Hereafter the convention Q = 10"(5„ will be adopted in c.g.s. 
units (e.g. L51 means luminosity in unit of 10^^ erg s""'^). In the initial state, the 
baryons are essentially at rest with respect to the central engine; 

• End of ejection (r = Aq): After the central engine constantly ejects energy 
(with luminosity L) for a time T, the fireball radius is r = Aq = cT = (3 x 
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101° cm) (T/l s); 

Coasting radius (r — Tc): As the fireball shell expands, the baryons will be accel- 
erated by radiation pressure. The fireball bulk Lorentz factor increases linearly 
with radius, until reaching the maximum Lorentz factor Fq, which is usually the 
dimensionlessentropy 77, but sometimes lower than 77 due to the limited radia- 
tion pressur^^EHl T he firebal l finally coasts with a constant Lorentz factor Fq 
at the coasting radiusPEHEH^ = TqAo ~ (lO^^ cm) (T/l s)(Fo/300) (this is 
the radius at which the entire mass of the fireball has achieved the coasting F). 
During acceleration, since all the materials essentially move with speed of light, 
the shell width in the fixed frame remains constant, A = Aq. The comoving 
shell width A' = FAq, on the other hand, increases with radius linearly. The 
above conclusions hold assuming the fireball shell evolves as a whole. However, if 
the mini-shells evolve separately, i.e., for well-separated GRB pulses, these mini- 
shells are likely to evolve independently, and coast at a smaller radius Vc = TghSo, 
where Tsh is the Lorentz factor of that particular mini-shell, and i5o = ct^ is the 
width of the mini-shells; 

Photospheric radius (r — rph)'. As the fireball shell expands, the photon number 
density and the typical photon energy drop. At a certain radius, the photons 
become optically thin to both pair production and to Compton scattering off the 
free electrons associated with baryons entrained in the fireball. At this radius, 
although much of the initial energy is converted to the kinetic energy of the 
shell, some energy will be radiated away as emission from rph with an approxi- 
mately blackbody spectrum (the fireball in the Big Bang has similarly a black- 
body spectrum, seen now as the cosmic microwave background emission from the 
last scattering surface). This is the first electromagnetic signal detectable from 
the fireball. This photosphere radius rph is usually above the coasting radius rc, 
with a temperature T = '?o('')9h/^c)~^^^, but could be below Vc if the initial fire- 
ball is clean enough (i.e. a large enough 77), in which case T ~ %. The typical 
value of the photosphere is Vpt ^ (lO^^ _ xq13) f^j. ^ ^ (^iqq _ j^qoO). A full 
discussion about various regimes of this "baryonic" photosphere is presented in 
Ref.Hni 

Internal shock radius (r = r^s): For an intermittent central engine with typical 
variability timescale of ty, the typical distance between adjacent mini-shells is 
usually also characterized as d = cty . Suppose that a rear shell moves faster than 
a leading shell, i.e., F(2) 3> F(i) ^ Fq, the (fixed frame) time at which the fast 
shell catches up with the slow shell is iis — d/(w(2) ~ W(i)) ~ 2T'^^^d/c, and the 
distance is ^ ksC ~ 2Tf^^d ~ 2Fgct^ ~ (6 x lO^^ ^m) (Fo/100)2(t„/0.1 s). 
At such distance s, mmi-she lls collide with each other and typically form strong 
"internal" shock J<^4 | 99 | l()U l 

Shell spreading radius (r = rj): For a shell with initial width Aq, assuming the 
Lorentz factor scatter AF within the shell is of order Fq, the shell starts to spread 
at a radius ~ FqAq, based on a similar argument as used to derive the internal 
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shock radiuJ^. The spreading radius of the mini-shells is roughly the internal 
shock radius. The spreading radius of the whole GRB shell is however much 
longer, typically r<; = (3 x 10^^ cm) (ro/100)^(r/10 s). Beyond the spreading 
radius Vg, the shell width starts to spread in the fixed frame, i.e., A ~ r/T"^ for 
r > Ts. 

• Deceleration radius [r — rdec)'- The fireball shell is eventually decelerated by 
the ambient medium (e.g. ISM). During the intial fireball- medium interaction, 
a reverse shock propagates into the fireball to stop it. Usually a deceleration 
radius (rdec) is defined as the radius where the reverse shock crosses the fireball 
shell. For a prompt fireball or a fireball whose duration is short enough, we have 
fdec = rr, i.e., the radius where the ISM mass collected by the fireball is equal 
to (I/Tq) of the firebaU rest mass, i.e., Mism ^ AM/Tq, where AM = MT is 
the total baryon loading of the fireball. For constant density medium, this radius 
is rr = imso/iTTnmpC^rl)^/^ ~ (2.6 x IQi^ cm) (i;..o,52/n)i/3(ro/300)-2/3. As 
long as the shell spreading radius is less than rr (or = fr/cTQ —> T), 
the fireball decelerates at rdec = rr- This occurs at an observer-frame time^^I 
tdec = tr ~ 5 (Ei.,n h?. /ny /HTn/3 00)'^^^(l + z) s. In literature this is usually 
termed the "thin shell" cas d I I Alternatively, if the shell is thick enough (e.g. 
for a long duration of fireball ejection so that T > tr), the decel eration radius 
moves further out to rr < rdec < rg- This is the "thick shell" cas J^^ I ^^ . As the 
fireball starts to decelerate, a strong external shock also forms and propagates 
into the medium. So the deceleration radius is essentially the initial external 
shock radius. 



Since the fireball moves essentially at the speed of light, the relevant fixed- frame 
(or central-engine frame) times corresponding to the fireball reaching various radii 
are simply derived by dividing the relevant distance with c. For example, the time 
for internal shocks to happen is about several hours while the deceleration time 
is about 10 days. If we imagine a GRB occurring in our neighborhood, say 10 pc 
from us, whose relativistic jet beam is exactly perpendicular to our line of sight, 
we would be able to follow the real time advance of the GRB emission and trace 
when the jet head reaches various radii. [For 9 — 71/2 and /i = 0, we have dt = dt 
in eq.l^]. However, when we see a GRB from cosmological distances, by definition, 
the jet is beamed towards us. The relativistic propagation effect (eq.P) squeezes 
all the timescales to within seconds. For example, the time delay between the the 
onset of the internal shock emission and the launch of the fireball (which may be 
due to the collapse event whose trigger time may be recorded by future gravita- 
tional wave detectors) is only tis ~ (ris/c)/2r§ ~ (18 ms) ris44(ro/300)"^(l -I- z), 
while the delay for the external shock emission is only tdec — tr — (rr/c)/2r§ ~ 
(5 s){E^so,52/nY^^{To/300y^/^{l + z) (for the thin shell case). In both expressions, 
the factor (1 -t- z) takes into account the cosmological time dilation effect. 



2, 2008 7:42 WSPC/Guidelines-IJMPA rev 



Gamma-ray bursts: progress, problems & prospects 27 

3.3.1. How common is the thick shell case? 
Although the thick shell case has been widely 

discussecPEnmnana, 

it is worth 

questioning whether a thick shell description is indeed relevant in reality. In all 
the current discussions, it is conventionally assumed that the fireball shell width 
is defined by the duration of the GRB itself, i.e. Aq = cT/{\ + z) (where the 
factor (1 + z) is again to correct the cosmological time dilation). By comparing 
Ts = TqAq wi th rr (or alte rnativelv comparing T with tr), one can define a critical 
Lorentz factoJnHClinni r, = 125(i;,«o,52/n)i/«(r/100 s)-3/8[(l + z)/2]3/8. For 
Fq > Fc, which is not difficult to satisfy, the fireball is in the thick shell regime. 
Hydrodynamically, the reverse shock becomes relativistic in the thick shell regime, 
while in the thin shell case, the reverse shock keeps Newtonian until reaching mildly 
relativitic at the deceleration radius. 

However, the above analysis is based on the assumption that the central engine 
has kept an essential constant luminosity through out T . In reality, the GRB energy 
injection is likely to be intermittent, sometimes with a broad gap between emission 
episodes. In such cases, it makes no more sense to take Aq = cT . Instead, one should 
separate the whole duration into several well-defined emission episodes, i.e., broad 
pulses, and treat these sub-shells to decelerate independently. One therefore has 
several consecutive thin shells, but rarely a thick shell. This situation has actually 
been indicated in GRB 990123, the only burst from which the rising lightcurve of 
the optical flash was caug In the standard afterglow model the optical flash 
is attributed to the reverse shock emission (see H5.1|l . The optical flash peak time 
(which corresponds to the time when the reverse shock cross the shell, ix = idee) is 
at 50 s after the burst trigger. The GRB duration, however, was 63 s, which exceeds 
fx- This already is incompatible with the thick shell regime, since in theory, ix 
should not be less than T . Furthermore, the rising lightcurve is very steep, which 
is consistent with a thin shell model, but inco nsistent with a thick shell model 
which predict a much shallower rising Inspecting the GRB lightcurve in 

detail, we find that the gamma-ray lightcurve of GRB 990123 is well represented 
as consisting of two components, a first more intense one which lasts less than 50 
s, and a second less intense one with even shorter duration. Thus, it is reasonable 
to regard 50 s as the deceleration time of the first shell in the thin shell regime. 
More generally, we suspect that in most cases when a GRB runs into an ISM, the 
thin shell case is very common. The thick shell description, may be more relevant 
in the case that a GRB runs into a pre-stellar wind or a constant dense medium. 
In such a case, the critical time for the thick shell case is much shorter due to the 
high density of the wind at the deceleration radiu a ■ 

3.4. Relativistic shocks 

The collision between the relativistic fireball and the ISM leads to a relativistic 
forward shock and possibly a relativistic reverse shock as well. The internal shocks 
are usually at least mildly relativistic due to the large Lorentz factor contrasts 



2, 2008 7:42 WSPC/Guidelines-IJMPA rev 



28 B. Zhang & P. Meszdros 

between the colliding shells. 

For a rela tiyistic shock, when the upstream matter is cold, the shock jump 
condition giveJ^^ZOl 



n-2 = (4721 + 3)ni ~ 472ini, 

62 = (721 - l)n2TOpC^ ~ 72lf^2n^-pC^ — 472inimpC^. (3) 



The subscripts "1" and "2" denote the unshocked (up stream) and shocked (down 
stream) materials, respectively; n is number density, e is internal energy, both 
measured in the comoving frames of the fluids, and 721 is the relative Lorentz 
factor between the fluids 2 and 1, while the Lorentz factor of the shock front itself 
is 7s = a/2721. 

In discussing the interaction between two fluids, two shocks are formed simulta- 
neously at the instant of contact, which propagate into the two fluids, respectively. 
There are four regions for the two fluids separated by the forward shock, the con- 
tact discontinuity, and the reverse shock. For the GRB external shock case, these 
four regions ar e (1 ) unshocked ISM, (2) shocked ISM, (3) shocked shell, and (4) 
unshocked shelt^. The shock jump condition (eq.|3]) also applies at the reverse 
shock. A final condition 62 = 63 (since p = e/3 for relativistic fluids, and equal 
pressure is required at the contact discontinuity) and the fact that the fluids move 
at the same speed across the contact discontinuity finally close the problem, and a 
simple analytical description about the system is 

availablfPI. 

Similar analyses can 

be applied to collisions between the mini-shells. Each collision is accompanied with 

a pair of internal shocks propagating into the two colliders. 

There can also be more complicated cases, e.g. involving three (or more) shell 

interactions. This happens in the phase when a fireball shell is already decelerated 

by the ISM, while a trailing shell (which could be ejected by the central engine 

at a later time, or be ejected by the engine at essentially a same time as the 

leading shell but with a much lower Lorentz factor) catches up with the decelerated 

leading shell and collides with it. This is a typical problem in the study of GRB 

afterglow evolution when there is additional energy injection (see t l3.7.3ll . and the 

941 

recent observed "step-like" optical afterglow lightcurve of GRB 03032& provides 
a strengthened motivation to study such a theoretical probleni. We have performed a 
detailed analysis of the three-shell-interaction hydrodynamic d l ^^'^ . In this case, there 
are altogether six regions separated by three shocks and two contact discontinuities. 
By applying shock jump conditions to all three shocks (and noticing that the leading 
shell is hot) , one can derive a self-consistent solution only when the relative Lorentz 
factor between the trailing and the leading shell exceeds a critical value defined by 
the energy ratio between the two shells. Otherwise, the injection is only mild, with 
the injected material connecting to the initial shell as a single long-lasting thick 
shell. 
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3.5. Synchrotron emission 

The GRB prompt emission is clearly non-thermal, and so is the afterglow emission. 
The most natural mechanism for non-thermal emission is synchrotron emission, 
i.e. emission from relativistic electrons gyrating in random magnetic fields. The 
question of whether the GRB prompt emission (e.g. the 7-rays) is definitely due 
to synchrotron is still subject to debate (see M.',^ . However, the synchrotron shock 
model is widely accepted as the major radiation mechanism in the external shock, 
which is thought to be responsible for the observed broad-band afterglows. 

There are three major assumptions that are adopted in almost all the current 
GRB afterglow models. Firstly, electrons are assumed to be "Fermi" accelerated 
at the relativistic shocks and to have a power-law distribution with a power-law 
index p upon acceleration, i.e. N{Ee)dEe oc E~PdEp Thi s is consistent with cur- 
rent shock acceleration numerical simulation j ' U45 | 14ij | Secondly, a fraction 
(generally taken to be <,1) of the total electrons associated with the ISM baryons 
are accelerated, and the total electron energy is a fraction eg of the total internal 
energy in the shocked region. Thirdly, the strength of the magnetic fields in the 
shocked region is unknown, but its energy density (B^/Stt) is assumed to be a frac- 
tion €b of the internal energy. These "micro-physics" parameters, p, (Ce) and e_B, 
reflect our inability of tackling the problem, whose values are usually fitted from 
the datal30 131 132 133 134 

There are several critical energies in the power-law distribution of the elec- 
trons. For p > 2 (which is consistent with numerical simulations, and seems to 
be consistent with most of the observational data), a lower limit is set by the re- 
quirement that the average energy density in the shock-heated region is j2in2'mpc'^ 
(eq.[3), which reads 7™ = g{p){mp/me){ee/^e)T ^ 310[5(p)/(l/6)](ee/Ce)r, where 
dip) — (p~ 2)/(p — 1) with p = 2.2 adopted, and we have simply redefined 721 as F. 
According to the standard synchrotron emission theory the radiation power of 
an electron is ~ {'^li)aTC'-i1{B'^ I^-k) cx 7^, so that high energy electrons "cool" 
more rapidly. For a continuous injection, as is the case in an afterglow (i.e. the for- 
ward shock keeps plowing into the ISM), there is a break in the electron spectrum 
at 7e = 7c, above which the electron energy spectrum is steepened due to cooling, 
i.e. N{Ee)dEe oc E^P^^dE^- This energy is time-dependent, which is defined by 
equating the comoving dynamical timescale of the blastwave {t' ~ Ft, eq.|2]) to the 
cooling timescale of the electron (t^ = 7emeC^/Pe)^^ Finally, the maximum en- 
ergy of the electrons (7m) is defined by equating the typical comoving acceleration 
timescale, t'^^^ ~ 27rri/c (where is the Larmor radius), with the shorter one of 
the dynamical timescale and the cooling scale. For electrons, the latter is relevant, 
which results in jm ~ {^e/arBy^^ ~ 5 x 10^(5/1 G)"^/^, where e is electron 
charge here. 

The typical observed emission frequency from an electron with (comoving) en- 
ergy 7emeC^ and with a bulk Lorentz factor F is = Tj'^{eB /2Trmec). Thus, three 
critical frequencies are defined by the three characteristic electron energies. These 



February 2, 2008 7:42 WSPC/Guidelines-IJMPA rev 



30 B. Zhang & P. Meszdros 




Fig. 11. The afterglow synchrotron spectra and lightcurves for the simplest fireball blastwave 
model (from Ref.86). 



are Vm (the injection frequency), Vc (the cooling frequency), and vm (the maximum 
synchrotron frequency). In the afterglow problem, there is one more frequency, Va, 
which is defined by synchrotron self-absorption at lower frequencies. So the final 
GRB afterglow synchrotron spectrum is a four-segment broken power law^^-^ sep- 
arated by the typical frequencies z/q, Vrn. , an d Depending on the order between 
Vm and z/c, there are two types of spectrJ^ (Fis. 111(1 . For < ^^c, which is called 
the "slow cooling case" , the spectrum is 



(z./z.™)-(f-l)/2 



Va<V <Vm 
Vc<V< Vm 



(4) 



For Vm > which is called the "fast cooling case", the spectrum is 



V <Va 
Va<V <Vc 

Vm<V <Vm 



(5) 



There are several more complicated regimes, involving self-absorptiorfSnHESI In the 
above expressions, the normalization factor is calculated by multiplying the t otal 
number of radiating electrons Anr^ni/S by the peak flux from a single electro 
which i s only the function of B and is independent of the energy (7e) of the 
electroiP™. 
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One direct consequence of sync hrotron emission is that the emission from an 
individual particle is polarized ^^ l Theore tical predicti ons about the polarization 
degree of afterglow emission have been ma Due to the random nature 

of the post-shock magnetic fields, the polarization is likely to be largely averaged 
out, and only a small degree of polarization is left. The most likely situation which 
can result in detectablehnear polarization in an afterglow is thought to be around 
the jet-break timJSZHESl jn which case a collimated jet and an off-axis line of 
sight conspire to produce an asymmetry which can lead to net polarization. In 
particular, for a l arge eno ugh offset of the line-of-sight, a 90° change of polarization 
angle is predictecl2l2E^ around the jet break time for the simplest uniform conical 
jet model. Current afterglow (optical) polarization observations detect low-degree 
(< 5%) of polarization emissiorr^^^ but the data are too sparse to adequately 
test the concrete models. 



3.6. Simplest afterglow model 



The simplest afterglow mode l i v l < " I v^ l g^ l vv l yO I is based on several assumptions, 
which minimize the complications. These assumptions include (a) isotropic fire- 
ball; (b) constant ambient density (ISM); (c) impulsive injection in the fireball; (d) 
relativistic fireball; (e) synchrotron emission of the electrons; and (f) constraints on 
the microphysics parameters (e.g. no evolution, p > 2, etc). Under such conditions, 
the fireball Lorentz factor evolves with radius r (or t) and with observer's time t as 

r cx oc t-^^\ r cx t^^'^ . (6) 

The first scaling is valid for an adiabatic evolution of the fireball in which the energy 
E oc nr^j'^ is constant, which is generally valid at late epochs (later than hours) in 
all afterglows, and is also valid at earlier epochs for many afterglows. The second 
scaling is valid for the evolution of a "radiative" fireball, iii which the to tal energy 
in the fireball decreases prominently due to radiation los J^^'^ l '^'*'^'^ ! ^'^^ ! ^'^ ^, while the 
momentum cx nr^j conserves. This extreme radiative regime is only valid for the 
fast-cooling regime when € p. ~ 1 ■ T-M^.A^^.t;- coii dition appears unlikely, according 
to the current afterglow fit J ^ ^^^^ 1 ^ ^ 1 ^ ^j^ l ^ ^j^? ! ^ ^j^ l ]\/[ore generally, a blastwave may 
not be strictly adiabatic or radiative. A reasonab le tre atment is to adopt a quasi- 
adiabatic evolution with small radiative correctiorPZ^I 

With a certain dynamics (e.g. adiabatic evolution), one can quantify the time 
evolutions of the critical frequencies, i^rm fc, i^a, as well as the normalization Fu^m, 
parameterized in terms of the burst properties (i.e. E^so^ n and luminosity distance 
di) and the micro-physics parameters (i.e. p, ^e, £e, and e^). Although the scalings 
on all these parameters are the same for different works, different treatments of 
normalization easily cause a factor of a few difference in the coefficients. Here we 
take the latest value of the coefficientJS, in which a typical value p = 2.2 has been 
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adopted. The dependence of unknown is also incorporated. 

lyrn = (6 X lO^s Hz) (1 + zy/'Eli\eJ^,)h]i\t/l day)-^/^ (7) 
i/c - (9 X 10^2 Hz) (1 + z)-i/2g^3/2^-i^-^i/2^^^^ day)-i/2 (g) 

z^a - (2 X 10^ Hz) (1 + z)-\eJ^,)-h]i\'/'E}f (9) 
F,,„, = (20 mJy) (1 + z)e]i\^/^ E^2dl%. (10) 

At a given time a comparison of Vm and Vc shows whether the flow is in the slow 
cooling or fast cooling regime, and using eqs.Q) and lO one can calculate afterglow 
lightcurves for a particular band (i.e. fixing a particular observation frequency). 
The different segments on the lightcurves have different temporal decay indices that 
correspond to different spectral regimeJSBI. The commonly used temporal indices 
include the following. In X-rays, shortly after the trigger, the temporal index is es- 
sentially constant, i.e. (2 — 3p)/4, which is steeper than —1 for typicalp values. This 
is consistent with the observations. In the optical band, the lightcurve first rises as 
oc i^/^, peaks around hours, and then decays with an index — 3(p — l)/4 ^ — 1. Al- 
though the rising lig htcurve has not been firmly detected (but ha; been inferrecina) 
due to the technical limitations, the falling part of the lightcurve is generally con- 
sistent with the observational data. The rising lightcurve may be buried by the 
contribution of the reverse shock emission at earlier times (as might have been the 
case of GRB 990123 and GRB 021211), but there is a large region of parameter 
space for which this part is not expected to be buried and should be observabl e-^-^ . 

A convenient way to test various model regimes is to perform simultaneous 
measurements of both the temporal and the spectral indices. Writing F^it^v) oc 
t^v^ , the relations of a and /? in various regimes are listed in Table 1. Also listed 
in the table are the jet, wind, and p < 2 models, discussed below. 

More sophisticated model takes into account the radiation of the whole bulk 
of the blast wave, which is usually modeled as a Blandford-McKe(jS2 self-similar 
profile, as well as the emission contributions from the so-called equal-a rrival-time 
surface, which is typically is egg-shaped or pear-shape d'^^^^l''^^'*^!''^^^!'"^ ^. These ef- 
fects result in smooth transitions around the spectral and temporal breaks, which 
provide a better model fit to the afterglow data. A complete description is presented 
in Ref . EHl 

Another important ingredient of the standard afterglow model is the emission 
component from the reverse external shock. During the early phase of the fireball 
shell - ISM interaction, a reverse shock propa gates into the fireball shell itself while 
the forward shock propagates into the ISm'''^'^. The reverse shock heats up the 
shell and accelerates electrons which emit synchrotron radiation as well. The shock 
is short-lived and ends when it crosses the shell. After this shock crossing time the 
electrons cool, leaving a rapidly decaying emission component. Due to the high par- 
ticle density in the shell compared to that in the external medium (e.g. the ISM), the 
peak frequency of the reverse shock is in the optical/IR band at the crossing tim e ty . 
Thus the reverse shock emission usually leads to an optical flasl t*^''^ ' ''^^''*^'*^ ' -^^-^^^ l -*--*-^ ! 
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Table 1. Temporal index « and spectral index /3 in various afterglow models, the convention Fi, oc f^v^ 
is adopted, from Refs. 85, 86, 110, 289, 326. The assumption Va < min(i/m,^'c) is made. (Under certain 
conditions, e.g. for the wind fast cooling case in some limited regime, the higher Ua case is relevant^*^, 
so that the values collected here are no longer valid). The jet model applies for the sideways expanding 
phase, which is valid for both ISM and wind cases and is usually in the slow cooling regime. 
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The early optical afterglow data of GRB 990123, GRB 021004 and G RB 021211 
are consistent with the reverse shock interpretatioj ^^'^ l ^^^ l ^^'^ l ^^^ l ^^" I ^Q » ^^^ '. We 
discuss in more detail the reverse shock and early afterglow signatures in ^5.11 

3.7. Additional features of realistic afterglow models 

The six assumptions made to define the simplest afterglow model are introduced 
mainly to simplify the problem. In reality, they are not necessarily satisfied, or not 
in all cases (although in many cases, they work surprisingly well). Violations of one 
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or more of these assumptions lead to second-generation or modified fireball models. 
Such effects are in fact naturally expected, and they can describe the physical prob- 
lem more realistically without the need for other ad-hoc assumptions. In this sense, 
they are also part of the standard fireball paradigm. The studies of the following 
effects comprise a majority of the recent GRB afterglow theoretical modeling effort. 
Because of space limitations, we will not go into the details of each of them, but 
rather refer the readers to the original papers that discuss these effects. 

3.7.1. Jets 

The "isotropic" assumption is one of the likeliest to be unrealistic, and investiga- 
tions of the consequences of assuming that GRBs are coUimated relativistic flows, 

i.e., jetJEHnna, 

have been fruitful. There are two quantities colloquially (and con- 
fusingly) referred to as the "beaming factor" in the GRB problem. One is the 
geometric beaming factor, which in the simplest jet model is just the opening angle 
of the jet, 6j. The second is the relativistic beaming factor, i.e., the emission of 
a particle or an object that moves with a Lorentz factor T is beamed into a cone 
with opening angle of l/F. More accurately, it is useful to refer to the first one as 
the collimation factor (or angle), and to the second one as the relativistic beaming 
factor. Jets give rise to an interesting interplay between the two effects. Initially 
the jet is ultra-relativistic, with l/F < 9j. An observer on the beam only receives 
information from within the relativistic light cone and has no knowledge about 
whether outside this cone the emitter is radiating or not. The description of the 
dynamical evolution is therefore equivalent to the isotropic case. As the jet slows 
down, eventually the relativistic beam becomes wider than the geometric beam or 
collimation angle, i.e. l/F > 9j. Two effects come into play. First is the edge effect, 
i.e., the observer starts to feel a deficit of energy per solid angle (which is the key 
parameter in afterglow theories, denoted as Eiso/An in the above notations). Sec- 
ond, the causally connected region starts to extend to the whole jet cone around 
the same time, and can keep expanding sideways. This means the jet can start to 
expand sideway^^. The times for the two effects to take effect are close to each 

othejnn, 

or may comcid depending on the assumption of the unknown ex- 
pansion speed. In the asymptotic regime for times much longer than the jet break 
time, the dynamics is 

F oc exp(-r/0 cx t^^/^ r (x t° (11) 

(where / = [Ej /{4TT/3)nmpC^]^/^ is the so-called Sedov length at which the collected 
ISM rest mass energy is equal to the jet energy Ej itself), so that the temporal 
dependence of various critical frequencies as well as the temporal decay indices 
all differ from the isotropic case (Table 1). In the asymptotic regime the post-jet- 
break optical lightcurve should have cx t^^ ~ t^'^, much steeper than the isotropic 
case (cx t^^)- Therefore a jet break is the natural interpretation for the apparent 
steepening observed in many GRB optical afterglows. 
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Although the behavior in the asymptotic regimes are well known, the detailed 
behavior near the jet break involves various complex effects including the jet spread- 
ing hydrodynamics. Many efforts have been made to model the smoothness of the 
111 11911^ 9^fi 

jet breaks iiozou^ j^^^ consensus is yet to be achieved even for the sim- 

plest uniform jet case. Hvdrodvnam ical num erical simulations are needed, but so 
far such studies are still preliminarj . Furthermore, GRB jets may not be 

simply uniform, but may be structured^^^"^^^. This could bring more interesting 
effects in modeling (see H4.4l for more discussions). 



3.7.2. Non-uniform external medium: winds and bumps 

The GRB ambient density may not be uniform. A well-discussed scenario involves 
an external medium produced by a stellar wind from the massive star progenitor, 
such as, e.g., a Wolf-Rayet star. The wind can be modeled to first approximation as 
ejecting material with a constant mass loss rate (M) and velocity (v^). Conservation 
of mass gives a p = Ar^^ density profilJ^, with A = M /AttVw = 5x 10^^^* g cm^^. 
Such a density profile changes the fireball dynamics into 

^ocr-l/2a^"l/^rcx^l/2, (12) 
so that the lightcurves for the homogeneous 

fireballHSEHEOEHi, 

for the co llimated 

jetPniMi, as well as for the early afterglow involving reverse shock emissio 
are all modified accordingly (Table 1). It is noticeable that most of the c urrent GRB 
afterglow data are consistent with a constant density external n^ediun^™ 
although a handful of bursts could be well modeled by the wind mode F^'^ l I '^' ■ ^ 1 
An external density j urnp w ould be expected at the interface of the stellar wind 
and the ISM outside of ilpSl which causes a distinct afterglow signature "'^''^ 
The ISM itself may also ha ve ji^fflj.^Y . §"^9^^ ^^^o^S- which would also add imprints 
on the afterglow lightcurve J'^^'^^ l '^'^^ l '^^1" ^^^. The consequence of a sudden drop in 
the density profile has also been discusseo^QS.. 



3.7.3. Post-injection and variable injection 

The fireball deceleration time depends on the ambient density, but the duration 
of the energy injection into the fireball is determined by the central engine behav- 
ior. These two timescales are independent of each other. Thus it is possible, or 
likely, that the central engine is still active while the afterglow starts. Thisgives 
rise to a post-injection of energy into the fireball, causing "refreshed shocks'*^, so 
called because the kinetic energy of the late-arriving material (catching up with 
the decelerated initial material) revitalizes the external shock. The injection could 
be either in a continuous form^''' '^^^^ '^^^ ^^■^ wi th modu lations superimposed, 
or discrete, e.g. through disrupted, collisional eventJ^^E^. The injection energy, 
both early and late, could be either in the form of a Poynting-flux-dominated flow 
or kinetic-energy-dominated shells ^ —'. However, the main characteristic of the post 
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injection signature is that the flux level is systemat ically increased after each injec- 
tion event without resuming the pre-injection level ^'l , a signature that has been 
seen in GRB 03032&3Q5.. 

Post-injection effect changes the energy per solid angle along the line of sight, 
and influences the fireball dynamics. If the fireball is not uniform in the angular 
sense but is patchy, the fireball deceleration will cause bright or dim spots to enter 
the relativistic beam, which equivalently changes the energy per solid angle of the 
fireball along the line of sight Both this model, and the varying density profile 
model, were proposed to interpret the recent observed optical lightcurve wiggles in 
GRB 02100 ^ '^' 1 '^^ 1 X-ray data could disentangle the two scenarios, since above 
the cooling frequency Vc emission is in the density-independent regime. For GRB 
021004, variability was also found in X-rayJ^, which lends supports to the patchy 
beam model. 



3.7.4. Relativistic to Newtonian transition 

A decelerating relativistic fireball eventually becomes non-relativistic (Newtonian) 
at the Sedov radius I = {3>E / AimmpC^Y^^ = (1.2 x 10^^ cm) (E'sa/n)!/^^ ^j^g 
collected ISM rest mass energy is eg ual t o the fireball enCT In the Newtonian 
phase, the fireball dynamics evolves a J^^ I ^^ I ^"'' l ^"^ l ^"^ l '^l" l 

?; (X cx r cx ^2/5. (13) 

The temporal decay indices ar — (15p — 21)/10 for Vm < v < Vc, and — (3p— 4)/2 
for > i^c, which for typical values of p is steeper than —1 as expected for 
the relativistic isotropic fireball, but is flatter than —2 as expected for the post- 
jet-break case. Therefore the evolution of a coUimated jet involves a steepening 
around the jet break (typically days to w eeks) and a later flattening as the fireball 
transmits into the Newtonian phasJ^^. This later transition time is In ^ l/c = 
(450 days)(i?52/n)^/^, which is longer than a year for typical parameters. Usually 
this is of observational interest in the radio banc . In order to interpret some shal- 
low lightcurve steepenings with relativistic-to-Newtonian transition effect (rather 
than the jet effect), a very dense medium (n '-^ 10^) is needed^. A self-consistent 
analytical description of the relativistic-Newtonian transition is proposed in Ref. 

ma 



3.7.5. High energy spectral components 

Besides synchrotron emission, there are other mechanisms giving rise to high energy 
spectral components which must be operative at some level, and which may have in- 
teresting observational consequences. A straightforward and widely-di^scussed com^^ 

ponent is the synchrotron self-inverse Compton (IC) componenll ^ l '^4 | 7-'^K^^^K^^^ I '^^ l '^^'^ l '^^4 | 315 | 31B | 317 | 318 | 319 | 18t] 

The IG effect plays two roles in studying GRB afterglows. First, electrons cool both 
via synchrotron and IC, so that IC potentially influences the value and evolution of 
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v^. Defining Y = Ljc/Lsyn, it is foundSHthat F = [-1 + (1 + ir]ee/eBy^^]/2 ~ 
{•qee/ esY^"^ for Y ^ \, where 77 is the overaU (synchrotron plus IC) radiation 
efficiency. The condition for IC cooling to be important is essentially rje,, > e^. 
The second role of IC is that it forms a second spectral component, extending 
beyond the high end of the synchrotron spectrum. To first order, it can also be 
approximated as a four-segment broken power law, separated by three critical fre- 
quencies v]p = Irn^a, i^jn — Im^^m, and I'l'^ — ^^Vc- In re ality the spectrum is 
more rounded and the power-law breaks are not as sharj^^l The high energy end 
of the spectrum is usually the Klein-Nishina limit . The flux normalization can 

be derived through Fl'^jFi.^rn ^ (16/3)crTCe?i»', which is < 10^^ for typical pa- 
il Sfil 

rameters and most times of interest) The condition for the IC component to 
stick out above the synchrotron component is defined by Fl^{i'l'^) > F^{i'l'^) for 
slow-cooling, and F^'^J v}^) > F^{i'lf;^) for fast-coohng, both relating to the same 
physical condition °" l Both this IC-emission-dominated condition and the above 
IC-cooling-dominated condition require cb ^ eg, although they are in principle dif- 
ferent from each other. It is worth noticing that current afterglow n iodeling suggests 
that the shock condition ^ eg is common among various burst J^ -*^ ! -*^'^'^ ! -*^^^ ! ^^^ !, 
so that IC is potentially important in afterglow physics. 

When observing in a particular fixed band, the IC component peak would sweep 
across the band at a later time than the sync hrotro n. This leads to a distinct bump 
signature in the hghtcurve. The IC peak time isP^t/c = (3.4 days) (ee/0.5)°-^^e^°^2 
£.e^^E^2°-°^n-°-^^{l + z)°-^'^i^:[s°-^^. The higher the energy band, the earher the bump 
shows up, and a denser medium helps to ease the bump condition needed to detect 
the IC component. A GeV a fterglow bump should be common for the currently fa- 
vored shock parameter J^^^ I In the X-ray band, an IC bump is expected to emerge 
on the power-law decaying lightcurve around a couple of days if the medium density 
is moderately dense (sa y, n > 5)317 319 186^ such a bump has already been 
detected in GRB 00092^123. 

Besides the IC component, several other high energy spectral components are 
also expected, which are related to hadronic processes. It is believed that the rel- 
ativistic shocks also accelerate protons besides accelerating electrons. The protons 
also radiate via the synchrotron mechanism, and can also interact with photons (e.g. 
synchrotron photons from the electrons) to produce pions and muonJ^2E22l222l 
The neutral pions decay into gamma-rays direct ly, a nd charged pions and muons 
also emit gamma-rays via synchrotron radiation I However, these components 
may not contribute significantly to the afterglow emission for the currently-favored 
shock parameters, and the parameter space regime for these components to be 



important is smalJ^^- 



3.7.6. Microphysics 



Although an electron distribution index p > 2 is consistent with numerical simula- 
tions of shock acceleratiorti^^lMSiMSl^ and is generally consistent with the observa- 
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tional in some bursts (e.g. GRB 010222) the temporal decaying slopes 

both befor e and afte r the "jet break" are too shallow, so that one seems to require 
1 < p < ^SSE23, In such a regime, the maximum energy power is distributed 
towards the high energy end of the electron spectrum, bringing some new features 
for afterglow emission. A possible mechanism for producing such distributions is 
given in Ref. IH87I and a study of this afterglow regime is presented in Ref. \IV2()\ the 
relevant conclusions being also collected in Table 1. 

Essentially all current afterglow model fits assumed non-evolution of all the 
shock parameters, p, e^, cb- In principle, these may change, but this is hard to 
quantify, and the effect may be degenerate with other effects so that one might 
never be able to disentangle them. Detailed afterglow fitting seems to be compatible 
with the model that one or more such parameters evolve with time . There 

might also be a gradient of magneti c fie lds behind the shock, but the data suggest 
that fields do not decay very rapidljE^. 



3.7.7. Pair formation, neutrons, grains and other effects 

There are a number of other physical processes that would modify the simplest 
afterglow model. 

First, the prompt gamma-ray form a radiation front which moves ahead of the 
blast wave and interacts with the ISM before the blastwave starts to decelerat e ' I 
The gamma-rays back-scattered by the medium can interact with outgoing gamma- 
rays and generate electron-positron pairs. The pairs enhance the opacity, and for 

a moderately dense external medium, the process leads to a run-away pair-loading 

'^'-iO '^^1 '^^9 '^^'-t 

process that modifies the blastwave dynamics considerablj 

Second, it is likely that there are free neutrons entrained in the fireball^ESaSSi. 

If neutrons and protons coast at essentially a same speed (which is not the case 

171 ^ '-i 7 

under certain conditions ' ) , the neutron shell would lead the proton shell 
when the latter is decelerated by the ISM. The comoving neutron decay time 
scale is ~ 900 s, so in the fixed frame, the distance where neutrons decay is 
Rp = CT^Tn — (0.8 X 10^^ cm) (r„/300). This leads to a neutron-decay trail ex- 
tending to a distance Rtraii ^ lORfj ^ 10^^ cm. The trailing proton blastwav e wo uld 
interact with this neutron trail and form a possible bump on the lightcurvJ^^. 

Third, if GRBs originate in star forming regions, the existen ce of dust grains 
would cause a new emission component, in the form of a "dust echo'^^^^Conversely, 
dust may be de stroy ed by the reverse-shock-induced UV- optical flaslPSi, or by the 
X-ray afterglow I The small- angle scattering of X-rays off the dust grains, on the 
other hand, may cause a soft X-ray bump on the afterglow lightcurvc^— . 

Finally, as supp orted by acc umulating evidence, at least some GRBs are associ- 
ated with 

supernovPEllSEni. 

Therefore, a supernova lightcurve, usually signaled 
by a reddened bump, is superposed on the optical afterglow lightcurve (as well as, 
in some bright cases, by a distinctive supernova spectrunJ^^l4% . 
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4. Problems 

A cartoon of a GRB fireball is shown in Fig. EH The theoretical model we have been 
discussing so far only includes those aspects that have been extensively tested, i.e., 
the afterglow part. This arises in the outer part of the fireball, from which most of 
the information is retrieved. This picture is widely accepted in the GRB community, 
and is used as a standard theoretical framework to confront with the afterglow data. 
In this section, we zoom in towards the GRB central engine, and discuss several 
aspects of the model on which consensus has not been quite reached. The closer 
towards the central engine, the less we know about the physical processes going on. 
Nonetheless, theorists have utilized the limited information available to construct 
toy models to tackle the problems. For some of the "problems" discussed below 
there exists a leading model to interpret the phenomenon (e.g. the internal shock 
model for GRB prompt emission). However, there is (are) other competing model(s) 
that is (are) still widely discussed, and there are continuing debates on some of these 
issues. 

A first major problem is the nature of the GRB prompt 7-ray emission itself. 
The origin of the gamma-rays has been debated since their discovery. Although the 
discovery of the afterglows settled the GRB distance issue, eliminated many GRB 
models, and achieved a consensus on attributing the afterglow radiation mainly 
to synchrotron radiation from a blast wave, there are still several variants within 
the current generic fireball scheme whi ch are argued about for interpreting the 
limited GRB prompt emission datd242J xhe major uncertainties include (a) the 
unknown fireball content (e.g. how important are magnetic fields in generating 
GRBs) and the unknown fireball energy dissipation mechanism (shocks or magnetic 
dissipation such as reconnection) , (b) the unknown location where the GRB prompt 
emission is produced (i.e. internal or external), and (c) the radiation mechanism 
(e.g. synchrotron radiation or other mechanisms such as Comptonization). These 
three aspects of the problem will be discussed below in H4. 114.31 We then move 
on to discuss the global GRB jet structures f ii4.4|l . progenitors fi i4.5|l and central 
engines ( H4.6(l . Finally we discuss GRB environments f H4.7|l and uncertainties in 
the shock physics ( H4.8|l . Some of these questions have also been discussed in Refs. 

[El El 



4.1. Fireball content: kinetic energy or magnetically dominated? 
4.1.1. Internal shock model 

The leading model for the GRB prompt 7-ray emission is the internal shock modei^, 
although there are various problems with it which are unresolved. The main mo- 
tivation for this model is to explain short time (down to ms) time variabilities^, 
more efficiently than e.g. external shocks^''' (c.f. Ref. I106|l . In this model, as dis- 
cussed in ^ it is assumed that the fireball behavior is essentially hydrodynamical, 
i.e. baryon or kinetic energy dominated, and magnetic fields are neglected when 
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considering the dynamics, and are only introduced when discussing synchrotron 
radiation through the parameterization of es- As discussed in HH.HI if the central 
engine ejects energy intermittently in the form of mini-shells, due to non-uniformity 
of the shell Lorentz factor, these mini-shells would collide at a typical distance of 
2TlcU (6 X 1013 cm) (ro/100)2(t„/0.1 s). The collisions produce a pair of 
shocks propagating into both shells, heating them and accelerating electrons (and 
protons). With some magnetic fields in the shell (either carried from the central 
engine or generated in situ, e.g. by turbulent dynamos or instabilities), these elec- 
trons emit synchrotron radiation in the gamma-ray band after Doppler boosting. 
This is taken to be responsible for the observed GRB emission. In order to have the 
gamma-rays from an optically thin medium (to avoid thermalization) , only internal 
shocks whose radii are above the baryonic photosphere {rph, see contribute 
to the observed emission. There is another relevant issue currently widely consid- 
ered. The gamma-rays generated from some closer-in internal shocks may have a 
large opacity to 77 interactions leading to pair production. The pairs generated at 
small radii form a Thomson scattering screen for the gamma-rays, and the radius 
where this Thomson scattering optical depth drops to unity defines another "pair 
photosphere", rpair, which is usually above Vph, limiting the minimum variability 
timescale of the lightcurves"'^^"^ '^^'^ The optical depth in this pair photosphere 
may be self-regulated to be around a moderate value of a few, so that variabilities 
shorter th an t he scale defined by rpair are only smoothed rather than completely 
washed ou 

The internal shock model has been extensively studied by 

various group^raraaimi. 

Its most successful aspect is its ability to model the complex temporal profiles of 
GRB prompt emission lightcurves. Numerical simulations indicate that the ob- 
served temporal behavior essentially reflects the temporal behavior of the central 
enginJP^. A caveat is that the central engine needs to be "bare" , or have a channel 
leading out to a lower density environment where optically thin shock radiation can 
be observed . In the currently favored "coUapsar" scenario f H4.5(l for long bursts, 
the central engine is located deep inside a collapsing star, whose envelo pe m ay act 
as an additional agent to regulate the variability of the relativistic flow^. In any 
case, as long as the relativistic jet emanates intermittently from the central engine, 
internal shocks are likely to develop (a different situation arises if the jet is strongly 
Poynting- dominated, as discussed next). The combination of the internal shock 
emission and the pair photosphere screen can reproduce the -5/3 slope and the 1 
Hz brea k in the power density spectra (FDS^^i "'^'-''^ which are deduced from 
the datPMlHI. 

An important characteristic of the int ernal sho ck model is that 
outflows with lower Lorentz factors have higher Ep 

vaiuejnma. 

This is contrary 

to most other models, and against the simple intuitive expectation that lower T 
(dirty) shells have a smaller Lorentz boost so the emission should be softer. The 
reason for this behavior in internal shocks is mainly that the low F shells col- 
lide at closer distances relative to the central engine, where the magnetic fields 
are stronger. A major reason why this model is widely studied is that there is a 
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clear operational theoretical framework (e.g. shock jump condition, equipartition 
parameters, synchrotron radiation, etc.) upon which model simulations could be 
performed (in contrast to the Poynting-dominated model discussed below). 

There are several problems or caveats about the internal shock model. First, 
the emission energy has to be extracted from the relative kinetic energy between 
the colli ding: shells, but the radiation efficiency is typically small, e.g. around 
In order to achieve a high GRB efficiency, as required by 
the fact that the pro mpt gamma-ray energy is or the same order of the afterglow 
kinetic enerr^^^^^ some novel suggest ions have been made, including non-linear 
dissipatioJS"! and quasi-elastic collisions'^^. In any case, a large relative Lorentz 
factor dispersion in the flow is required2diS4£J.0SI, A recent study indicates that 
by taking into account neutron-dec ay in the fireball, the efficiency of the internal 
shock mech anism is further reduced ^^ 'I Second, the BATSE bright burst spectral 
sampl^^3 suggests that, at least within the same burst, the peak photon energy 
Ep distribution is narrow. This is hard to achieve within the internal sho ck rnodel 
unless one invokes a strong bimodal distribution of the Lorentz factorj ^^ l ^^*^ ! 
which lacks straightforward physical origins. Although none of the above criticisms 
has displaced the internal shock model from its position as the leading paradigm, 
these are important issues to consider. For this reason, other models are worth 
investigating as well, some of which are discussed below. 



4.1.2. Fireball with strong magnetic fields 

There are several motivations for considering strong magnetic fields in the GRB 
problem, as opposed to the weaker fields needed to produce the observed radiation. 
(1) Electromagnetic energy is "clean" in the baryon load sense, and can propagate 
in vacuum. Since a GRB fireball typically requires a very small baryon loading 
to achieve the high Lorentz factor needed to solve the compactness problem, a 
Poynting flux potentially make s it easy to transport a large amount energy with- 
out carrying much baryonJ^^ (2) Current GRB central engine models commonly 
invoke a rapidly rotating black hole circulated by a debris torus, or a millisecond 
neutron star. During collapse, magnetic flux conservation naturally gives a field 
of order 10^^ — 10^^ G, as is observed in radio pulsars. With rapid rotation near 
the bre akup frequency, magnetic fields are likely to be magnified via an a — f2 
dynamo '^ 1 to achieve ^ 10^^ G or higher. Such high fields are also argued to be 
possible for direct collapse of a magnetized white dwan=S^. A strongly-magnetized 
rapidly-rotating central engine is a likely scenario; (3) Magnetic fields are a pos- 
sible agent to tap energy from the two energy reservoirs in the engine, i.e., the 
gravitational energy of the torus, and the rotational energy of the black hole or the 
neutron stai-'^^-'^ ^'^'^ ^'^'^ 156^ Another energy extraction mechanism, i.e., the neu- 
trino pair production process {ly'D — > e+e~)^, is found to just barely power GRBs 
when a beaming correction is taken into accounlPSSE^Iso that the magnetic power 
is at least helpful in meeting the GRB energetic need g ^ (4) Magnetic fields are 
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helpful in coUimating Other arguments in favor of the magnetic mech- 

anism include its ability to alleviate the inefficiency problem of the int erna l shock 
modePSSI^^ 

and the possibility of achieving narrow Ep distributionJ^^I, There 
are also three observational facts that support a strongly magnetized central en- 
gine (although the flow is not necessarily completely Poynting flux dominated). (1) 
The strong gamma-ray polarizatioiJ33 (cf. Refs. El il may indicate a str ongly 
magnetized central engine, either in pure Poynting-flux-dominated forrrP^, or in 
conventional hydrodynamical form but with a globally organized magnetic field 
configuration'^'-'^ '^"-''^; (2) Modeling the reverse shock emission of GRB 990123 indi- 
cates that the reverse shock region should a nchor a s tronger field (by a factor of 15 
in strength) than in the forward shock regioJEDEnni. A similar conclusion may also 
apply to GRB 02121 1-'^'^'-' These may indicate that the fireball contains a large 
portion of "primordial" magnetic fields carried from the central engine, although 
the flow is not completely Poynting-flux dominated (otherwise there should be no 
strong reverse shock'^^-'^). (3) The non-detection of a bright photospheric thermal 
component may need a significant fraction of energy being stored in the magnetic 
form, so that the photosphere tempera ture could be reduced (i.e. by a factor of 
(1 + cr)-!/-*) to evade detectioJ^^S^ (but see ReflMJ. 

The evolution of a fireball with large magnetic content, and how such a fire- 
ball i^ ^y giy^ rise to prom pt gamma-ray emission has been reviewed by various 
author J ^^^ l ^^^^ l '^^^ l '^^^ l '^^ , but so far there is not a standard framework, as is the 
case for non-magnetic fireballs. The magnetic fireball evolution is now an inter- 
play among three components, an electromagnetic component, an internal energy 
component, and a kinetic (bulk) energy component. On the other hand, traditional 
fireballs only invoke the latter two, with the internal component dominating in 
the beginning, converted into the kinetic form via radiation-pressure-driven accel- 
eration, and partially re-converted back to the internal form via shock dissipation. 
The total amount of energy of these two components is essentially conserved except 
for the radiation losses. To simplify the magnetic case, one can regard the sum of 
these two components (internal and kinetic) as one single component, through a 
parameter a which is the ratio between the energ y den sities in the electromagnetic 
component and in the internal-|-kinetic componen ff^ I At the beginning of the fire- 
ball evolution when the kinetic component is negligible, a is just the ratio of the 
cold component (in the form of Poynting fiux) and the hot component (in the form 
of photon-pair fireball generated via neutrino annihilation)*^. During the coasting 
regime, when radiation-driven acceleration is completed, the ratio a is simply the 
Poynting-flux- to- kin etic-energy ratio, as has been widely discussed in pulsar wind 
nebula theorieJ^Sl In the GRB problem that we are interested in, a is likely to 
decrease with radius, since the rapid spin of the central engine greatly eases the 
conditions for magnetic reconnection and other instabilities to develop, so that the 



'^We note that the coexistence of both components is natural when a cataclysmic collapse event 
forms a rapidly rotating, strongly-magnetized compact object. 
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magnetic field energy is dissipated". About a half of the dissipated magnetic en- 
ergy is converted to internal energy, and the ot her half energy is used to accelerate 
the fireball through magnetic pressure gradient! '^^^^ I '^'^'^ l '^'^^ l '^^^ l If the dissipation 
radius is smaller than the photosphere radius, the internal energy is also converted 
to the kinetic energy eventually. However, if the dissipation radius is above the 
photosphere radius, the internal energy would be radiated event ually v ia the non- 
thermal electrons accelerated during the reconnection even y" ' '' ' ' l '" ' 'v l vvv' l This is the 
main source to interpret the GRB prompt emission in this model. 

The inclusion of electromagnetic fields greatly complicates the GRB problem. 
For not very high a flows, one simplified treatment is magnetohydrodynamical 
(MHD). The adoption of such an approximation is justified when the electromag- 
netic fields are fully coupled to the fluid, i.e., no dissipation is allowed. As discussed 
above, a GRB magnetic fireball needs to be intrinsically dissipative in order to in- 
terpret bursts, thus casting the MHD approach in question. Nonetheless, the MHD 
method is still used to solve the first-order problcm"^^"^ 354^ assuming that the ra- 
diation loss is only minor, and that a is not very high. For a high-enough a (e.g. 
o'>,o'c2 ~ 200), the MHD condition breaks down globally at a distance r^jjD = (1.8x 
10^^ cm) Lgf [cr(l-|-cr)]-i/2(t„/l ms)roj ~ (1.8 x lO^^ cm) Llf^^a^^ty/l ms)roj, 
which is defined by the c ondition that the real plasma density is lo wer than what 
is required for corotatior|S3^3 i.e. the Goldreich- Julian densitjE^. Here ac2 
is a critical value (at the relevant radius, here it is t'mhd) above which r^^-^-^ < 
rdec, i.e., the MHD condition globally breaks down before the fireball starts to 
decelerati^SI Although the condition for this to happen is stringent'^'''' 342^ when 
it happens there might be a global magnetic field dissipation region within the 
radius range r^^^^ < r < Vdec, in which electrons are randomly accelerated by 
turbulent electromagnetic fields, giving rise to high energy radiatio n (al though the 
typical energy may not coincide with the observed GRB emission^ '-^ ' "^ l In reality, 
such a high a value may not persist out to such a large radius, especially when 
reconncction-induced dissipation presumably take place from the very beginning of 
the flow'^^'' "^'^'^ 'i70 Thus the MHD description may be generally valid to reveal 
the first order phvsicJ^ '^^^ l '^' ^"^ ! Howe ver, the formalism is still rather complicated 
even in the MHD formalisrrPSniSHI When the fluid inertia is negligible, as is the 
case in the high-cr limit, an even simpler treatment, i.e., the force- free approxima- 
tion, could be introduced^^^l^^. Using this simple formalism, the ev olutio n of a 
magnetic bubble is delineated, and is found (surprisingly) to be similaJ^^to the 
self-similar evolution of the hydrodynamical blast waves- ^. An interesting finding 
is that a point explosion gives rise to structured angular distributions of energy and 
Lorentz factoi'^^''', i.e., e{d) oc (sin^)^^ and F cx (sin6')^^, which nicely matches 



'^In pulsar wind nebula theories, on the other hand, it has been a long-standing problem to solve 
the so-called ir-problem, i.e., there is no obvious reason to reduce a from a pulsar wind (slow 
pulsars compared with GRB central engine) from a very high value (10'*) to a low enough value 
(< 1) required in the reverse shocks as inferred from the observations. 
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one version i 



of the structured jet (see more discussions in i j4.4(l . __ 

An important characteristic of a high-tr flow is that strong shocks cannot develofP^. 
This is because the shock frame pressure is dominated by the magnetic field pres- 
sure, and not much internal energy is available for radiation. The condition for 
strong shocks is cr < 0. j '^^^ l As a consequence, the internal shock mechanism has to 
be replaced by something else if the GRB outflow is st rongly m agnetized (e.g. a > 
1). Besides reconnection- induced magnetic dissipatiorPSaESHI other me chani sms 
include Comptonization of the photospheric emission via Alfven turbulencJ^^ and 
synchro-Compton radiation in a large amplitude electromagnetic wave as the mag- 
netic fireball is decelerated by the ambient ISlvl^ZS in all these models, however, 
the fireball is assumed to be converted back to hydrodynamical in the deceleration 
phase i n order to match with the achievements of the hydrodynamical afterglow 
theorieJSMESZI. 

Such a transition, although seemingly ad hoc at first sight, is in 
fact reasonable since such low-cr flows have been observed in pulsar wind nebula, 
even though the initial pulsar wind is deflnitely a high-cr flow. 

The development of magnetic-field-powered GRB models is still preliminary. 
The difficulties that hamper its progress lie in the intrinsic complication introduced 
by electromagnetic fields so that it lacks an operational scenario to work on. The 
situation has been changing recently "^"-"^ '^'^^ "^"^^ '^'-''^ '^'-'^ S67 q^^^ potential problem 
for the magnetic model is that an extremely clean fireball involves much fewer 
baryon-associated electrons. Giving the same amount of dissipation energy , the 
typical energy of emission tends to be much harder than the sub-MeV band ^^ I 
unless secondary pairs are involved in the problem, which o n the ot her hand tend 
to smear out sharp variabilities observed in GRB 

lightcurveJ^32E3. 



4.2. GRB location: internal or external? 

A second uncertainty in understanding GRB prompt emission lies in its location 
within the fireball. The non-thermal spectra suggest that the location should be 
above the photosphere radius Vph, below which the emission is opaque. Afterglow 
radiation limits the prompt emission radius to be smaller than the d eceleration ra- 
dius rdec- Generally there are three suggestions in the literaturJ^^l (x) the exter- 
nal modek^^ ^^'^ ^^"^ ^^"^ '^'^^ that suggest GRB prompt emi ssion occurs upon 
strong deceleration of the fireball at rdec', (2) the internal model j^ '^ l '^'^ l ^^^^^ l '^^^^ l '^^'^ l 
that suggest GR B occurring at a radius Vph < r < rdec', and (3) the (inner- 
most) modelJ^^^''^2JI3 342]iti2at suggest GRB occurring right above the photosphere 
{r>rph)- The current consensus is that the third emission component at most con- 
tributes to the final emission only partially, while a debate between the first two 
components (external or internal) is still going on. 

The main focus of the debate is the GRB variability. For those bursts con- 
taining chaotic, spiky pulses, an internal scenario seems more natural to interpret 
the d ata. The variability could be due to the intermittent behavior of the central 
enginJ^, including that induced from black hole accretion instability, the intrinsic 
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instabilities involved in the jet propagation within the stellar envelopJ^, and the 
intrinsi c cha otic behavior of magnetic reconnections and instabilities for the high-cr 
scenaricP^. However, whether an external model can be ruled out is still uncer- 
tain. For a burst running into a homogeneous ISM, if the central engine lifetime 
(teng) IS shortcr than the dec eler ation time tdec (thin shell), the observed duration 
of the burst should be simp tang ~ ^dec/r^c, which is causcd by the effect of 
angular spreading f H3.2(l ^^. In such a case, an external shock model only produces 
one single smooth pulse. The caveat is that the external medium could be clumpy, 
and if the scale of the clumps d is smaller than r^ec/r, variable lightcurves could 
be produced. An argument against this scenario is that the process is likely to be 
inefficienJSS the efficiency being estimated as 77 ^ ty/T <C 1 to first-order. However, 
taking into account the angle-dependent flux effect -'^^'^ 376^ ^j^g observed variability 
may still be reproduced with a higher efficiency. So solely from this fact, the two 
scenarios can not be differentiated. The attractive feature of the internal scenario 
is that the variability arises more naturally, while the external shock model has to 
introduce an additional assumption, i.e., a (very) clumpy medium. It would then 
be worth checking the further consequences of such a clumpy medium. Evid ence for 
non-uniform medium has been suggested for GRBs with X-ray line feature^^^Sl 
and for GRB 021004 which s hows a wiggling afterglow lightcurve^^^ and multi- 
component absorption featureJl^. However, in most afterglows, a clumpy medium 
is not required, and limits may be posed on the dumpiness. It is worth modeling 
the prompt GRB emission as well as the afterglow emission within the external 
shock model in a unified manner, and to make use of any statistical correlations 
between the GRB variability and afterglow irregularity, so as to validate or falsify 
the external shock GRB model. 

Other clues can also help to differentiate the two scenario s. (1) The correlation 
between the "waiting time" and the amplitude of the next puls o^^ I favors an engine- 
dominated (internal shock) model; (2) In the external shock model, one expects 
pulse widths spreading with tim C^^^ , which is not apparent from the data. However, 
this objection may be circumvented by taking into account the noisJ^^^ (3) The 
absolute value of F enters the problem in the external shock model, while only the 
relative value of F enters the problem in the internal shock model. As a result, 
some correlations are expected in the external shock model but not necessarily in 
the internal shock model. For example, a br ightness-spectral hardness correlation 
is expected in the external shock modePl, which is consistent with the data. A 
brightness-duration anti-correlation is also expected, but this is not apparent from 
the data. A statistical study about various correlations among duration, flux and Ep 
in the external shock modeP^ suggests that the model predictions are compatible 
with the data, but the results do not necessarily rule out the internal shock model. 
For the internal shock model, the burst duration is expected to be independent 
with flux and spectral behavior. 

Eventually, close monitoring of gamma-ray and X-ray emission at the transi- 
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tion time of prompt emission and early afterglow would be able to clearly identify 
whether or not there is a distinct afterglow component to stick in at a later time. 
This would eventually settle down the internal-external debate. This is in principle 
doable in the S wift e ra,. Preliminary evidence of two distinct components is available 
for some bursts"^. 

One point (which is sometimes overlooked) is that the internal and external 
shock scenarios are not mutually exclusive. Both are part of a generic cosmological 
firebah*^^ '^'^'^ model. It is likely that both components coexist in all fireballs, 
with different components dominating in different bursts. For example, in some 
parameter regimes, the internal shock radii are beyond the external shock (deceler- 
ation) radius, so in this case, no internal shock is expected. There is a_small fraction 
of BATSE bursts that show single-peak FRED-like tempo ral p rofiles^ . These events 



poralpr 



are most naturally interpreted as of external shock origii 

In the magnetically-dominated fireballs, there are also an internal emission 
component'^^'^ ^'''^ and an external emission componentf2^^«2£L. A highly variable 
lightcurve is possible for both cases. 



4.3. GRB emission mechanism: synchrotron, or others? 

The leading radiation mechanism introduced to interpret GRB prompt emission 
is synchrotron emissiorP"^ ^'^ '^^ This is the most natural mechanism, which 

is found successful to interpret afterglows and many other astrophysical phenom- 
ena. Many o f the observ ed GRB spectra are found to be consistent with this 
interpretatioiPSOEHIEHS. 

However, the synchrotron model encounters difficulties 
in explaining several observatio na.1 fa cts . The most notable one is sometimes called 
the synchrotron "line-of-death'^^^Sni { ^^ the low energy photon number spectral 
index should not exceed -2/3 (or the spectral index should not exceed 1/3). 
A good portion of GRBs seem to violate this limilJ^^ in the BATSE database. 
Proposals to solve this apparent inconsistency include critical discussions of the 
sensitivity of BATSE to giv e acc urate slopes at such low energieJ^^, introducing 
small- ang le "jitter" radiatio rPSl, anisotropic electron pitch angle distribution and 
synchrotron self-absorption^^5I2I2l^ as well as photospheric emissiorx"^^ . A nothe r ar- 
gument raised against the synchrotron mechanism is the cooling problerrP^. The 
synchrotron cooling time scale is very short, so that the prompt emission should be 
in the "fast-cooling" regime with very low cooling frequency, and an Fi, spectral in- 
dex of -1/2 (or photon number spectral index -3/2) should obtain, either for the low 
or the high energy spectral indices, but this is not apparent in the data. This crit- 
icism may be circumvented by taking into account more complicated acceleration 
and cooling processes . 

Although the synchrotron model continues to be the main paradigm, some other 
alternative radiation mechanisms have been proposed to interpret GRB prompt 
emission. The most natural extension is to involve Compton scattering. If the 
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typical synchrotron emission frequency is well below the sub-MeV band^, then 
the GRB prompt emission could be due to synchrotron self-inverse- Compton (IC) 
emissiorl^lSl. Caveats about this suggestion include, (a) a strong IC component 
indicates a higher radiation energy density than the magnetic field energy den- 
sity, so that a positive feedback results in increasingly stronger higher-order IC 
components before the Klein- Ni^hina cutoff, and consequently greatly increases the 
total energy demand of GRbJS^, (b) a strong IC component requires es/ee <C 
2^317 319 186^ while in internal shock scenarios, eg can not be too small given a 
strongly-magnetized central enginJ^^l (c) an IC mechanism tends to furt her in- 
crease the Ep dispersion due to the high power dependence of Ep on l'^^, which 
within the internal shock synchrotron model is already problematic without intro- 
ducing bimodal distributions of F't/"'^ o^o 

Alternatively, models involving Comptonization of thermal or quasi-thermal 
particles have been proposed. A saturated Compton cooling model seems to be 
able to fit the observed GRB spectra weli^^^*^^. It is unclear how the model pa- 
rameters needed for the fitting could be generated naturally in the fireball model, 
although attempts to transplant this mechanism to the internal shock scheme has 
been madJ^^. Comptonization off a very dense photon bath as the GRB mecha- 
nism was also proposcG , with the soft photon field being provided by the trapped 
photons within the massive envelope funnel of the progenitoJ^^. A proton-induced 
pair-photon synchrotron cascade model, in analogy to a nuclear pile reaction, was 
proposed recently^S^j,^ which can reproduce a narrow Ep distribution around 0.5 
MeV (electron rest mass). For the usual luminosity and bulk Lorentz factors, and 
when cosmol ogica l redshift correction is incorporated, this seems to be consistent 
with the datEpSl. 

The rep ort of a measurment of strong linear polari zation in gamma-ray prompt 
emissior@2l, though subject to debate and confirmatior03^3 potentially holds im- 
portant clues for understanding the GRB radiation mechanism. Models pro ducing 
high gamma-ray polarization include Poynting-flux dominated model a I , mod- 

i nyolying hydrodynamical shells with entrained globally structured magnetic 
fieldJSSHEMl^ and models involving off-beam observers for a narrow jeJSlElil. 
Models involving inverse Compton scattering with offset beaming angles can also 
give rise to large degrees of polarization"^^*^ 397 .. 399 .ii ^^9 8^ Additional information is 
needed to differentiate these possibilities. 

Within the standard synchrotron model, there is a s traig htforward way to an- 
alyze the Ep distribution within various fireball models'^ In general, one can 
write 

Ep - T-il{neB/mec){\ + z)-^ - (2 x lO"* cV) {TB)-fl{l + z)-\ (14) 



'in fact, in order to raise the synchrotron typical energy to the BATSE band, some parameters 
need to be pushed to extremes. For example, eg needs to be close to unity, and (the f raction of 
electrons that are shock accelerated from the shocks) needs to be sometimes less than unitjEHZj^HI 
so that the energy per electron is raised. 
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where F is the bulk Lorentz factor of the firebaU, 7e is the typical comoving Lorentz 
factor of the electrons, B is the comoving magnetic field, and z is the redshift. In 
principle, B could have two origins, one carried from the central engine, and an- 
other generated locally (e.g. in the shocks), although recent evid ence sug gests that 



the engine-origin magnetic field may be the dominant componenlP^ELl Both field 
components, once formed, have the same dependence on distance from the central 
engine, i.e., the transverse part goes a,s B (x r~^, if both a and es do not change 
substantially with radius . The r— dependence of 7e essentially determines the 
Ep distributions for the models. For the internal shock model, since 7e is determined 
by the relative Lorentz factor between the two colliding shells, it should not have a 
strong dependence on r, so tha.t Ep cx {TB){l+z)-^ cx L!/'^J~}(} +z)~^ ■ The depen- 
dence on r reflects a strong dependence of F (Ep cx jS2EIISl which is the intrin- 
sic reason why the internal shock model cannot reproduce a very narrow E' y^ distri - 
bution within a same burst (unless bimodal distribution of F is introduced I * ) . 
For a magnetic-dominated acceleration model, on the other hand, 7e should only 
depend on the local field B, generally as 7e oc B^^/^ if Lamor acceleration is in- 
volved. This nicely cancels out all the r-dependences, so that Ep oc F(l + z)^^. 
Only the dispersion of F and z contribute to the final dispersion of Ep. This is 
one of the strengths of the magnetic-dominant internal models"^'^^. For the external 
models in both the low- and high-cr regimes, the case is more complicated with the 
interplay of the ambient density. The Ep distribution is found to be wide if indepen- 
dent distributions of luminosity and other parameters are assumec^-*^. A narrower 
Ep dist ribut ion is possible if some intrinsic parameter correlations are taken into 
accouni^^. 

Different dependences of Ep on various parameters for different models provide 
a framework on which the models may be tested with future data (Table 2)'^^'3. 
The parameters on which Ep depends in various models are in principle measurable 
with current and future observational facilities. Two major unknown parameters 
include the redshift and the bulk Lorentz factor. The former is measurable in large 
num bers based on expected localizations by the next generation GRB mission, 
scheduled to be launched in 2004. The latter could be al so m easured or 
at least constrained in the Swift era using the early afterglow dati^"^ (see more 
in ij5.1|l . or could be measured using the high energ y cut off of GRB spcctra^'-'2l 
when the next generation gamma-ray mission, GLAST^^, is launched in 2006. So 
there are good prospects for pinning down the location and mechanism of the GRB 
prompt emission as well as the content of GRB fireball through statistical analyses 
of future data. 



4.4. GRB jet: uniform or quasi-universal? 

The geometrical configurations are an essential ingredient in characterizing and 
understanding astrophysical phenomena (e.g. pulsars, AGNs, etc.). Evidence sug- 
gests that GRBs are very likely collimated. Understanding the degree of collimation 
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Table 2. The dependences of the GRB spectral break energy {Ep) on various parame- 
ters in the synchrotron emission model of various fireball variants''''^. Model parameters: 
L: initial total luminosity (not just the radiated luminosity) of the fireball; E: initial 
total energy of the fireball; F: bulk Lorentz factor at the radius of GRB radiation; ty : 
typical variability time scale; n: ISM density; z: redshift. 
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as well as the possible structure of the collimated flow would be very essential to 
understand the burst mechanisms and the true event rates. 

The main evidence for GRB fireball coUimation is provided by the achromatic 
steepening breaks in some GRB optical afterglow lightcurves. The simplest model, 
i.e., a conical jet with a uniform energy distribution within the cone and sharp 
energy depletion at the jet edgt-'^'^^ '^'^-'^ (Fig. [T^ . and how this model inter- 
prets the lightcurve steepenings, have been introduced in H3.7.1I In this model, the 
time for the break to occur roughly corresponds to a measure of the jet opening 
angle, if the observer's line-of-sight is in the jet cone and not too close to the jet 
edge. Data show that such inferred jet opening angles {6j) show a large dispersion 
among different bursts, and so do the corresponding isotropic gamma-ray emission 
energies (E^^iso)- However, both dispersions conspire in such a way that E^^isoO'j 
is essentially a constanl l'-^" ! "^" l Within the uniform jet model, the total energy re- 
leased in a GRB event can be estimated as — (i?^.iso/47r)27r(l — cosOj) x 2 — 
(1 — cos9j)E^^iso, where 27r(l — cos6'j) is the solid angle of a conical jet with an 
opening angle 6j, and the factor 2 takes into account the consideration that the jet 
is likely bipolar. At small angles, /& = (1 — cos^j) ~ sin^ 0^/2 ~ Therefore the 

above fact suggests a constant energy reservoir among long GRBs"^*^ '^^^^ One puz- 
zle posed by this is why a standard energy is collimated to rather different degrees 
among different bursts. 

An alternative, and in principle more elegant, model is to postulate a quasi- 
universal jet structure, i.e., there is a non-uniform distribution of energy per solid 
angle within the jet, so that all burst outflows could more or less retain a similar 
geometry, the diversity of observed jet break times and isotropic-cquivalcnt energies 
being caused by a different line-of-sight relative to the jet axis-*^^"^ ^'^^ ^'^'^ . If such 
a model can be constructed, GRBs might be considered to have not only a stan- 
dard energy reservoir, but also a standard geometric configuration. A lesson learned 
from AGN studies is that much of the apparent diversity in AGNs is simply caused 
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Fig. 12. Comparison between the quasi- universal structured jet model and the uniform jet model. 

by viewing-angle effects. Current unified AGN models try to find a paradigm in- 
terpreting various phenomena as being due to a standard configuration whicli is 
viewed at different angles'^'-'^ "^^^ . A unified picture for the GRB phenomenology 
would be similarly appealing^. A straightforward speculation is to introduce a jet 
configuration such that e{9) cx 6*"^, as suggested by the Ej,iso9j ^ const empir- 
ical law. Here e{9) is the energy per solid angle along the direction defined by 9 
(the angle between the viewing direction and the jet axis), which by definition, is 
equivalent to Mtt. Su ch a configuration has been discussed earlier to model 

afterglow lightcurveiS^EHZI but the viewing direction was placed close to the jet 
axis. By placing the viewing direction at arbitrary angles -'^'^'^ 136 ^ jt is found that an 
achromatic lightcurve steepening is naturally reproduced. However, the time for the 
steepening to occur now corresponds to the epoch when the line-of-sight Lorentz 
factor r{9v) is decelerated to a value below (rather than below l/9j for the 
uniform jet model, see Fig. El for a comparison between both jet models). 

^Before the standard energy suggestioiE^ it had been speculatetBSHl that GRBs may have both 
a standard energy reservoir and a possible standard configuration based on a handful of afterglow 
data. The geometric configuration suggested in that work includes three distinct components, i.e., 
two uniform cones (one narrow and one wide) and one quasi-isotropic component. There was no 
discussion about whether such a configuration is consistent with the afterglow lightcurves, and it 
seems that such a structure is likely to cause some distinct lightcurve signatures that may violate 
most of the data. It is , howe ver, likely to be consistent with the recent two-component jets as 
inferred in GRB 03032d4M2( 
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The e{d) cx 6^^ model (hereafter k = 2 power law structure model) is directly 
motivated by the E^^iso^] ^ const empirical law. From the point of view of re- 
producing jet steepening with the viewing-angle effect, more general types of jet 
structure, e.g. power law models with more general indices e{9) oc or even 
non-power-law structure, such as a Gaussian profile, may do the same job '^v l This 
is because well before the lightcurve "jet break" time, what is relevant for the dy- 
namical evolution as viewed by the observer is only the average energy per solid 
angle within the ^ 1/T{9y) cone, which essentially remains the same as the rela- 
tivistic beaming cone gets larger (given that [r(0i,)]^^ ^ 6y is satisfied initially). 
After the jet axis enters the field of view, eventually all the initial jet structure is 
expected to be smeared out due to energy redistribution and sideways expansion. 
So in both asymptotic regimes, a structured jet has the same temporal evolutions 
as the uniform jet model l Different jet structures only manifest themselves dif- 
ferently around the jet break time. The lightcurve of a fc = 2 power-law structured 
jet has been modeled by various authors -'^'^^ ^'^^ ^'^^ ^-'^'^ ^-'^-'^ ^21. It is found that 
the sideways expansion effect is not prominent, so that a model with t his effect 
neglected can still roughly reproduce the basic feature of the lightcurv^S^EHll 
For more general structure functions, the lightcurves are not adequately mod eled 
due to the complicated physics involved. Modeling Gaussian-type jet evolutior0^ 
revealed a substantial energy structure redistribution during the evolution. The 
resultant lightcurves are consistent with the jet data. 

A structured jet model has several clear predictions to confront with the ob- 
servations. (1) Because different luminosities in this model are caused by different 
viewing angles, whose probability is well defined, the structured jet model a spe- 
cific prediction about the GRB luminosity function for bursts in the same redshift 
power-law model, one has N{L)dL oc L'^-^/^dL, while for the 
Gaussian model, one has N{L)dL cx L~^d^^^. When the parameters are allowed 
to have some dispersion (quasi-universal picture), Monte Carlo simulations show 
that the above scalings still hold, with turnovers nea r the low- or high-luminosity 
ends caused by the variations of the parameterM In contrast, the uniform jet 
model has no predictive power concerning the luminosity function. Currently, al- 
though the sample of bursts with spectroscopic redshifts is too small to allow a direct 
luminos ity function study, th e luminosity functions derived using various statistical 
methodJ "^'1^^01^^1117^2221 

are not inconsistent with the model predictio ns, espe- 
cially for the quasi-universal Gaussian jet model with constant total energjE^I^^. 
(2) By taking into account the cosmological effects, the total number of bursts de- 
tected (regardless of redshif ts) a s a function of jet angle also has a clear prediction 
in the structured jet model I For the k — 2 power-law model, a distribution 
peak is predicted around 0.12 rad, which is in rough agreement with the current 
data. The peak is ex pecte d to move to larger angles when the detection sensitivity 
threshold is increase d I Again the uniform jet model has no prediction power on 
this. In the Swift era, both of the above predictions may be tested. 
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Besides the preliminary support suggested by the above two tests, the structured 
jet model may also be indirectly supported by other theoretical and experimental 
results. (1) Progenitor and central engine studies naturally give rise to jet structures. 
Simulations within the collapsar model naturally predict a jet structure after the 
jet penetrates thro ugh the stellar envelopJ^. Studies of magnetic-threaded ce ntral 
engine 

modelJSSl 

as well as the evolution of a Poynting-flux dominated flo-wl^^ 
result in jet structures in a natural way; (2) The interpretation of some GRB prompt 
emission empirical relations, such as the spe ctral lag - lun iinositv correlation , are 
consistent with the structured jet 

hypothesiJSHESninni^ 

so that a coherent picture 

is achievable for both the prompt emission and the afterglow emission within this 
theoretical framework. 

Several criticisms and caveats for the structured jet model have been raised. (1) 
The uniform jet model is very simple, and one can question the need to introduce 
jet structures unless it is necessarj 'tUI. On the other hand, a non-uniform jet struc- 
ture is a natural expectation in any jet model which arises from realistic physical 
processes. Even if the structured jet model so far does only the same job as the uni- 
form jet model, it is still worth exploring, and should be explored, since it is more 
physical and has more predictive power; (2) The k = 2 po wer-law model predicts 
some anomalous signatures inconsistent with the dat a I However, this may be 
caused by the assumption of an unphysical singular point in the jet struct ure, and 
is removable when more realisti c Jet structure (e.g. Gaussian) is considere d ^ (3) 
Recent optical polarization dat a I indicate a near 90° change of the polarization 
angle , whi ch is inconsistent with the model prediction in the k = 2 pow er-law jet 
model I However, the data accumulated are so far not conclusiv e ! ° l 

An important contraint on the jet structure is posed by the recent discovery that 
an Ep (X {EisoY^"^ correlation earlier proposed im GRB app ears now to extend all 
the way down to energies char acteristi c of XRF j ^ / M g l This relation, together 
with the Eiso oc correlatiorP^ESEI immediately leads to the inference Ep oc OJ^ , 
which poses important constramts on both the simple k = 2 power law model and 
the on-beam uniform jet mode^^. This is because Ep varies by two orders of mag- 
nitude from GRBs to XRFs, so that even if XRFs correspond to events viewed at 
the equator (or isotropic events), GRBs would have to be events viewed within the 
1° viewing angle (or corresponding to jets with opening a ngle sroa ller than 1°). For 
the fc = 2 model, this tends to greatly over-generate XRfJ^^EZI^ which are fo und 
to contribute ~ 1/3 of the total GRB/XRF population in the HETE-2 datJSHI. Por 
the uniform jet model, the very n arrow GRB jets contradict or do not address some 
important afterglow jet break dat j^'^^ l All current GRB/XRF prompt emission and 
afterglow d ata a re however consistent with a quasi-universal Gaussian-like struc- 
tured mode^^, with an angular structure e[9) ~ eq exp(— 6'^/26'o). In this model, 
GRB/XRF jets still retain a characteristic angle 6*0, with a mild structure inside and 
a rapid exponential decay outside. The XRFs are only those events with viewing 
angles 6'^ ~ (3 — 4)610, which greatly decreases the number of XRFs. Statistically, 
the jet parameters (e.g. the typical angle, total energy in the jet, etc.) are allowed to 
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have some scattej ' ■ ^ l Monte Carlo simulations indicate that with reasonably small 
scatter, a quasi- universal jet model not only solves the GRB/XRF population prob- 
lem, but can also reproduce the Eiso 9J^ corr elatioiMl^. More complicated 

jet structures are in principle possible, at least for some eventJiHuninni. These 

models along with the simplest Gaussian model should be more extensively con- 
fronted with the various data, including luminosity function, redshift distribution, 
opening angle distribution, etc! '' I I 

To conclude, evidence in support of long GRBs having a quasi-universal struc- 
tured jet configuration is mounting, but issues and contradictions remain. This 
paradigm will be more fully tested with the advent of the extensive data sets ex- 
pected in the Swift era. 



4.4.1. Other jet models 

Besides the uniform conical and the structured jet models, there are some other 
GRB geometric suggestions discussed in the literature. _____ 
One suggestion is that GRB jets might be 

cylindricaP2aS2ni, 

as observed in 

some AGNs. Current modeling assumes that the line-of-sight is on the jet beam, 
and no apparent jet break is predicted from the model, so that the model can not 
interpret the whole GRB phenomenology, but may account for a sub-category of 
the bursts. A cylindrical jet is howev er consistent with the MHD description of a 
Poynting-flux-dominated 

Another suggestion is that GRBs arise not fro m fluid flreball jets, but are rather 
"cannon balls" ejected from the central 

enginJlMEa. 

The picture is taken from 

analogy with some observations in micro-quasars and supernova remnants, although 
it is unclear how cannon balls could survive as compact entities against instabilities 
if they are accelerated over a large dynamical range of radii and Lorentz factors. The 
model has several distinct predictions (e.g. superluminal motion of the source, non- 
resolvable source image, and so on) that are different from those from the standard 
fireball jet models, which may be easily tested with the current and future data. The 
assoc iation of GRBs with supernovae has been quoted as a strong support for this 
nrodePM, but this piece of evidence is not exclusive, since GRB-SN associations 
are widely expected within the standard model as we L^^'^^'^ ' '^'^^ ' ^^^ ' ^^ ' ^*^ '. 



4.4.2. Orphan afterglows 

An important phenomenon associated with jets is the issue of "orphan afterglows" , 
i.e., the implication of low energy (X-ray, optical, radio) decaying transients which 
are not be detectable in 7-rays. The idea initially comes from the uniform jet 
modcx-'^'-'^. The starting point is that the initial GRB jet is sufRciently misaligned 
respect to the line-of-sight so that the observer misses the bright gamma-ray emis- 
sion. Later as the jet decelerates, it expands sideways and enters the line-of-sight. 
The observer is expected to see a steeply rising afterglow lightcurve followed by a 
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normal decaying lig The ev ent rates of orphan afterglows in vari- 

ous bands have been estimatecP^HESnESOl b^t so far there is no firm detection of 
orphan afterglows of any kind. 

The issue of orphan afterglows is complicated by two other considerations. First, 
there might be cosmological "dirty fireballs' or "failed GRBs" whose Lorentz 
factors are too small to allow transparent gamma-rays to be detected. These fire- 
balls, when decelerated by ISM, can also give rise to orphan afterglows, whose 
signatures are not easy to differentiate from those of the jet orphan afterglow^^. 
Some possible differences between these two types of transient events have recently 
been proposed I Second, if GRB jets are structure d I I, the orphan after- 
glows of the first kind would be greatly reduced since there is no sharp jet edge 
cutoff. The second-type (dirty fireball) may exist due to higher baryon loadings 
at larger angles from the jet axis. So far there is no detailed discussion about the 
orphan afterglow rates for the structured jet model. 

4.4.3. Nature of X-ray flashes 

So far there have been essentially four types of interpretations for the XRF phe- 
nomenon (see il2.1.1| in the literaturJ333, i.e. dirty fireballs or failed GRBs4M|3Z|432l, 
high redshift GRBijSZI, 

fireballs dominated by photosphere emissiorl^S 370 ^43' ^^^ 



geometry-related events, such as offbeam GRBs in the uniform jet modei^^^'^^ 
and large viewing angle events in the structu red j et modc^-'^"^''' The last model 
is likely in an anisotropic supernova explosion '^^ l The first two interpretations at- 
tribute the soft-faint nature of XRFs to the dispersion of two parameters, i.e. F 
and z, respectively''. In reality, Ep is determined by the combination of many pa- 
rameters (' H4.3(l . so that XRFs are hkely to be caused by t he di spersion of at least 
several of th em, as indi cated by Monte Garlo simulationJ^SI The photospheric 
interpretatio n * v i v f U l tf require the Ep distribution of the combined sample 

of GRBs/XRFs to have a sep arate component in the low energy regime. Recent af- 
terglow observations for XRfJ^^ pose important constraints on some of the above 
suggestions (e.g. high-z bursts and dirty fireball). In the Swift era, a large data set 
of both GRBs and XRFs, with photometric redshift measurements, should become 
available, which will eventually pin down the nature of XRFs. At present, the geo- 
metric model, espec ially the one involving structured jets, appears to us the most 
promising possibilit; 



:iaiiy 



4.5. Long burst progenitor: collapsar or supranova? 

Accumulating evidence suggests that long GRBs are associated with deaths of 
massive stars. The first cosmological progenitor scenario, i.e., the NS-NS merger 



''Note that the dirty fireball interpretation is inconsistent with the internal shock model, which 
predicts a higher Ep for a lower F (' i|4.c{i . It is however consistent with the external models and 
with the internal model involving a Poynting-dominated flowl242l 
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Fig. 13. Comparison of the geometric configurations of the collapsar and supranova models (from 
Ref.5). 

modems, is now disfavored for lone; although it is still a leading 

contender to account for the short, hard bursts (cf. Re f. 44 03. A variety of GRB 
progenitor models have been discussed in the literatur J-^^^ l and the require- 
ment for accretion onto a central object occurring in a neutrino-dominated regime 
significantly constr ains s ome of the scenarios"'^'^^. Although the formation ofGRBs 
from binary systemJ^^ or from gravitational energy loss in neutron star J'^^^ J have 
also been discussed, the leading scenario for causing a long GRB involve the core 
collapse of a single massive stai^^ "^^^ ^'^^ ^^"^ ^^'^ . A supernova (SN) 

explosion is naturally expected to be associated with the GRB (although in its 
first incarnatiorP^ this was thought to be a "failed" supernova, i.e. a core collapse 
which did not to achieve ejection of its stellar envelope). The launching of a GRB 
jet is widely believed to involve a black hole - torus system resulting from this 
core collapse, hence the GRB is either produced simultaneously with the SN, if 
the star collapses to a black hole promptly, or the GRB is delayed with respect 
to the SN (which would give rise to a "supramassive" , rapidly rotating neutron 
star as a first step, the neutron star later collapsing to a black hole in the second 
step, after loosing the angular momentum required to sustain the mass). Either a 
one-step or a two-step collapse are possible, depending on the mass and angular 
momentum of the progenitor core as well as the details of the collapse. If the delay 
of the second collaps e is not very long, e.g. from minutes to hours as suggested by 
numerical simulationJ^, the two-step collapse scenario is not very different from 
the one-step collapse, and both of these cases can be referred to as "collapsars'^^, 
e.g. of type I and II. On the other hand, an alternative suggestion is that the delay 
is long, e.g. from days to weeks or even months, and in this case the burst environ- 
ment is very different from the collapsar case, due to the role of the well-separated 
supernova shell and a central pulsar wina^*^ This latter case is referred to as 
a "supranova" . Figure 1131 is a cartoon picture for the geometric configurations in 
both models. 

The collapsar model has many merits in reproducing the data. First and fore- 
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most, it naturally involves a GRB-SN association, and predicts that GRBs are 
associated with star forming regions. The existence of the stellar envelope helps 
to collimate a je t with angular structure and can help to regulate the jet flow 
intermittently^. The duration of a burst is set by the fall-back timescale rather 
than the accretion timescalJ^, which is natural to interpret the durations of long 
GRBs. The requirements for the progenitor include a narrow range of specific angu- 
lar momentum (3x 10^^ cm^ s~^ < j < 2x 10^^ cm^ s~^) and poor metallicity. Such 
a constraint may lead to a rough standard GRB energy and an event rateroughly 
consistent with the true GRB event rate after coUimation correctionJ^^. All of 
these are consistent with the main-stream observational and theoretical progress in 
the GRB field. The interaction of such a jet with the stellar envelope should also 
make unique signatures in both electromagnetic and non-electromagnetic forms. For 
example, prompt X-ray and gamma-ray signals are expected when the i et breaks 



through the envelope, leading to precursors to the main buj-gijMIM^ 

no 

leading to short, hard bursts . The breakout of the jet cocoon is also a candi- 
date to interpret the X-ray emission line features hours after the burst trigger as 
have been seen in several bursts. Internal shocks within the jet before penetrating 
through the en velop e would lead to high energy protons to produce prompt TeV 
neutrino signal J I A collapsar may also harbor a long -lived central engine that 
continuously pumps energy into the GRB firebal l I I Recently, more realistic 
MHD simulations for the collapsar model have commence which will reveal 
the possible magnetic nature of the fireball. The difficulties of the collapsar model 
include the following. First, the required high angular momenta of the progenitors 
are difficult to achievJ^. Second, although the top candidates for coUapsars are 
those stars that have undergone intense mass loss before collapse (e.g. Wolf-Rayet 
stars), the expected wind environme nt for G RB afterglows^Q^S^ is not commonly 
identified in current afterglow studieJSHinil. 

The supranova model was initially introduced to alleviate the baryon-loading 
prob lem , since weeks to months after the SN explosion, the environment is relatively 
cleaiPl. Two other incentives were added later. One was its promise in interpreting 
X-ray line emission (and absorption) features. The conjectured supernova shell (or 
some torus-like or funnel-like remnant) located at 10^^ — 10^^ cm from the central 
engine (which was ejected by the progenitor days or months before, moving with 
a speed of ^ 0.1c) provides a large mass of heavy elements for the inferred pho- 
toionization (by the continuum emission of the GRB or the early afterglow) and 
recombination processes to produce a strong Fe lint"'^^'^ 4oO J23 which are thought 

to be responsible for many of the discrete X-ray spectral features detected in several 
burst J1«7 | 34 | 3,1 

(More recently, lower-Z elements such as Ca, N, S etc have also 
been reported in burstJSHE^ Such an interpretation of X-ray lines has some dis- 
tinct energetical advantage d ^ 1 but it requires high dumpiness in the shell. This 
type of model is not the only one able to explain the lines. A competitive class 
of models attributes t he sp ectral features to continuous photoionization by a long- 
lived central engine jell I or by a jet cocoon in the collapsar scenario, interacting 
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with the outer layers of the stellar progenitoi'^^l These models make the line at 
smaller distances (r ~ 10"'^^ — 10^'^ cm, hence they are referred to as "nearby" mod- 
els) and in a much higher density medium, require a small er iro n mass. Detailed 
modeling indicates that they are able to reproduce the dat M Other X-ray line 
models propos ed in clude irradiation of a pair screen caused b y the early afterglow 
back-scatterinJ3^ and the Cerenkov line emission mechanisnPSSI 

A second advantage of the supranova suggestion is related to the role of the 
central pulsar after the SN explosion; this can form a pulsar wind bubble which 
can modify the ambient medium from a pre-stellar wind to a medium with con- 
stant (lower) number density, consistent with the obser vations. H igh equipartition 
parameters (e,., e^) as inferred from afterglow modelin^^^^22l are also naturally 
interpretea^^^ 

The supranova model suffers several criticisms. First, it ta kes fine -tuning; to make 
a GRB months to years after the formation of the neutron staJ ^ ^^ I ^ . Second, lacking 
an envelope, the collapse may not produce a long burslP^ and th e cql limation 
mechanism is not clear. Third, detailed radiative transfer calculation^^^ indicate 
that a supranova model would reproduce the reported large equivalent widths only 
for very shallow incidence angles of the ionizing continuum, which is naturally 
expected in nearby models (involving ionizing the walls of a stellar funnel), but is 
less natural for an ionizing continuum incident on a distant supernova shell, where 
normal incidence is expected. Finally, the most severe objection to the supranova 
model is the fact that the association between GRB 030329 and SN 2003dh allows 
at most a very short delay (less than two days) between the SN and the GRB 

Aft 

events, and is compatible with both events being simultaneous . This rules out 
the supranova model at least for this burst, and supports the collapsar model. 
Unfortunately, no high spectral-resolution X-ray observation were obtained in a 
timely manner, thus missing an excellent chance to test whether a pre-ejected SN 
shell is a requisite to generate X-ray lines. 

Although the GRB 030329/SN 2003dh association greatly changes the balance 
between the collapsar/supranova debate, the issue of whether there is still a sub- 
category of bursts originating from supranovae (e.g. perhaps those which have 
strong X-ray lines) is still unsettled. Future observations can help to distinguish 
the two scenarios. First, both scenarios predict different electr omagnet ic signals. 
The collapsar scenario predicts X-ray and gamma-ray precursor^S^I^, while the 
interaction between the fireball and the SN shell in the supra nova model tends to 
produce enhanced, delayed GeV emission signal^^^EnnEnOl Second, the supra- 
nova model predicts more emission components and stronger fiux levels of high 
energy neutrino emission due to the existence of the SN shell which p rovides more 
target photons and protons for and pp 

interactioJSEMlIMr 

This leads to 

stronger neutrino signals both for individual sources and for the diffuse background. 
Detection/non-detection of such an excess emission may be used to prove/disprove 
the supranova model. Finally, suggestions to directly search for SNs that occurred 
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weeks to years before GRBs have been proposecP^, which may pose direct con- 
straints on the models. 

4.6. Central engine: what is behind? 

The most widely discussed GRB central engine invokes a central black hole and 
a surrounding torus. There are three ultimate energy sources, the gravitational 
binding energy of the torus, the spin energy of the black hole and the magnetic en- 
ergy. Two ways of extracting the accretion e nergy and black h ole spin energy have 
been c onsid ered, i.e. the neutrino mechanisni I I I I I, and the Blandford- 
Znaj mechanism. The former mechanism typically powers the conventional 

"ho t" fireball, while the latter invokes strong magnetic fields threading the black 
hol g I I I ' Iwhich typically would power a (at least initially) "cold" Poynting 
flux dominated flow. The identification of the content of the fireball f i|4.1|l and the 
mechanism of GRB prompt emission ( ij4.3|) would shed light on the mechanism that 
powers the central engine. 

Another type of GRB central engine often discussed in the literature is a mil- 
lisecond magneta J^i'^ ''^ '''^ 468.469 . 470 . 471.443.365 . 472. rj,^^^^ rapidly rotat- 
ing neutron stars with surface magnetic fields of order ~ 10^^ G and higher. 
These objects would be able to satisfy the essential conditions for a GRB cen- 
tral engine, e.g., energy, duration, variability, baryon loading, birth rate, etc, es- 
pecially when considering the possibility that initially there could be some tempo- 
rary toroidal fields with much stronger strength (10^^ Q^|4£lil44sML^ xhe ultimate 
energy sources include the spin energy of the pulsar and the magnetic energy. 
For a pulsar engine, there could be in principle two energy components, an initial 
prompt component (either via neutrino mechanism or via destroying the temporary 
toroidal field) powering the prompt GRB, and another long-term component due 
to spindown of the millisecond pulsar. The latter component continuously injects 
energy into the fireball through Poynting flux'^'^'^ '^^'L^^^ j or in the form of electron- 
positron pairs analogous to the pulsar wind bubble J^^ leaving well-defined bump 
signatures in the afterglow lightcurve J'^^^''^ l '^^^'^ l "^^^ l ^^'^^ . It also naturally provides a 
long-lived central engine, which may be the agent to continuously photo-ionize the 
Fe in the stellar envelope, giving rise to the observed X-ray emission features In 
the two-step supranova models, a pulsar is invoke d in the first step, which powers 
a pulsar wind to drive 

a magnetic-enriched bubb« 
A more exotic central engine mechani srn in volves a phase transition from normal 
neutron matter to strange quark matten^^. The process is likely in a detonative 
mode, leavin g beh ind a star completely composed of strange quark matter, called 
a strange staJ I The possibility of neutron star - strange star phase transitions 
powering GRBs has been discussed within the context of one-step^'''^ '^'''^ '^'^ or 
two-step processe^^^EZni A special category of such models lead to intermittent 
energy injecti on d ue to unstable photon decay, which is arguably a viable GRB 
central engin^SOl A common caveat about the strange star mechanism is that 
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there is no evidence yet about the existence of strange quark matter, while the 
existence of black holes and pulsars have been widely tested. 

4.7. Environment: what is in front? 

On the galactic scale, GRB host galaxies are broadly similar to the normal, star- 
forming faint field galaxies at comparable redshifts and m agnitudes . The distri- 
bution of GRB-host offsets from the galactic centeJ^Sl as well as heavy element 
abundance studieJi^ are fully consistent with a progenitor population associated 
with sites of massive star formation. 

On the parsec scale, issues about the GRB immediate environment include 
whether the medium is (quasi-) uniform, wind-like or clumpy, and how dense the 
medium is on average. Insights about these questions can be obtained through af- 
terglow l ightc urve modeling -'^'^-'^ ^'^'^ ^'^'^ -'^^^ and time-dependent absorption feature 
modelin g^ I, but the situation is still controversial. As mentione d earlier, the con- 
stant density medium model is consistent with most afterglow dat^J-^ ^'^^ l ^'^'^^ ^^, 
but the p oc wind model also works in some bursts "^^"^ "^^^ ii jg puzzling 
how different GRBs could have quite different immediate environments if they come 
from the same type of progenitor. There are so far no studies on correlations of 
the inferred GRB environment with other properties of the gamma-ray prompt and 
afterglow emission. The inferred mediun i densitv also varie s significantly among dif- 
ferent modelers, even for the same burslP^ ^'^^ l ^'^ l ^'^'^ l "^'^^'^ !, although there seems to 
be a trend towards favoring a universal m odera te-dense medium with n ^ 10 cm^'^ 
as being consistent with (most of) the datcPS. 

Suggestions that GRBs are embed- 
ded in dense molecular clouds have been mad^SHl which is consistent with the high 
hydrogen column density (Nh) inferred from some X-ray afterglow studie^^^Sl] 
A high-density afterglow model, on the other hand, although inferred from some af- 
terglow fits^^ may not be compatible with broad-band data-'^'^-'^ ^"^'^ or with the 
inferred standard energy budget for all GRbJ^. The dust grains expected in molec- 
ular clouds may reflect and irradiate the afterglow emission a nd for m a bright dust 
echo'^'^^, or a distinct IR signature detectable for nearby burst j^ ^. Du st extinction 
is also invoked as the cause of at least some optically dark burstJ^^. So f ar direct 
evidence of dusts is not yet collected, and the possibility of dust destructiorPSHSni 
further complicates the issue. In some bursts (e.g. GRB 021004), multiple spectral 
absorption systems and density bumps are identified, which refer to a rather non- 
unif orm, bum py environment from close to the burst to much farther away from the 
^^j,gjl59 |487|^ It has been sugge sted that time-dependent absorption features both 



in optica l "" I and in X-ray d lare powerful tools to study GRB environments, but 
the present data are too sparse to allow firm conclusions to be drawn. 

4.8. Shock parameters: universal or unpredictable? 

It is almost certain that afterglows are produced by coUisionless relativistic shocks 
energized by the GRB. The physics involved in these shocks is however poorly 
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known. A widely discussed scenario is that particles are accelerated via repeatedly 
crossing a shock front, a nd achiev e a power-law energy distribution through the well- 
known Fermi mechanisrrPSSEnOI TJ^e latest numerical simulations indicate that the 
resultant power-law index is universaf^^^ i,e., p ~ (2.2 — 2.3). This is in sharp 
contrast with what is inferred from afterglow fit data. With the simplest jet model, 
broadband modeling indi cates tha t p is quite unpredictable, ranging from ^ 1.4 to ~ 
2.8 among different 

burstjnnna. 

An important caveat for the modeling is that the 
present model only uses the simplest assumptions, e.g. uniform jets, non-evolving 
equipartition parameters, etc. The electron power-law index p is usually derived 
from the temporal decay index after the lightcurve break or steepening attributed 
to jet properties (i.e. F^, oc t~P). Whether the inconsistency between broad-band 
modeling and shock acceleration simulations is caused by the simplified GRB model 
or by the simplified shock acceleration model is not known, and developments in 
both directions are needed. For example, proposals for generating a flat electron 
index have been outlineo , while more complicated afterglow models (e.g . pos t 
injection, or refreshed shocks) could change the inferred p value considerabl y I 

Unlike the magnetic fields invoked in the magnetically dominated scenarios for 
GRB emission, which may be carried out from the central engine, the origin of 
the magnetic fields in the internal or external shocks needs to be addressed. The 
inferred fields are typically of strength ~ lO^'^ — 1. Loosely, these have been 
attributed to a turbulent dynamo mechanism behind the s hocks . More specific 
suggestions include two-stream instability of re latiy istic plasm^^ or the presence 
of an intrinsically magnetized ambient mediurrPSSI 



5. Prospects 

We have summarized the achievements (33 E^nd fJ^J and uncertainties ([Q in our 
current understanding of the nature of GRBs. Guided by past experience, we believe 
that our knowledge of GRBs will be further advanced in the coming years, driven 
by future observational breakthroughs. In this section, we attempt to foresee some 
of the possible milestones in the upcoming new epoch of GRB study led by NASA's 
two future missions. Swift and GLAST, as well as by experimental developments 
in new channels of GRB study. We will discuss how our knowledge of GRBs is 
likely to be extended in the temporal domain (i.e. early afterglows and the prompt 
emission - afterglow bridge, ti5.1|) . in the spectral domain (i.e. GeV-TeV emission, 
ii5.3|) . and into other, non-electromagnetic regimes, including high energy neutrinos 
f ^5.5l) and gravitational waves f tl5.6|) . We will also discuss how future developments 
could unveil the mystery of short GRBs f ti5.2|) and how GRB studies may become 
a unique tool to investigate the early universe f H5.7|) . Some of the topics covered in 
this section are also discussed in Refs. ^1 
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5.1. Early afterglow and reverse shock emission 

The discovery of GRB afterg has greatly extended our knowledge of 

GRBs, both in the temporal domain (hours to even years after the burst) and 
in the spectral domain (below the gamma-ray band, from X-rays down to radio). 
However, due to the long time scale of the mission alerts up to the ground, the 
slowness in slewing instruments and/or bad weather in ground-based instruments, 
in most cases afterglow observations have started hours after the burst trigger. At 
this phase, the afterglow blastwave has been decelerated and entered a self-similar 
regime, and the behavior is essentially determined by the total energy per solid angle 
in the fireball and the properties of the ambient medium. Some precious information 
characterizing the fireball, such as the initial Lorentz factor of the fireball and the 
magnetic content of the fireball, has at this stage been lost. In order to retrieve such 
information, very early afterglows need to be studied, which contain information 
on the emission from the reverse shock region during and shortly after the reverse 
shock crosses the fireball shell at the very beginning of the shell-medium interaction. 
As of October 2003, only four bursts (GRB 990123, GRB 021004, GRB 021211 and 
GRB 030418) were caught within less than 10 minutes after the triggerJ^^, thanks 
greatly to the growing number of robotic optical telescopes spreading all over the 
world. The situation will change soon following the launch of Swift in June 2004, 
which will automatically record essentially all the triggered GRB early afterglows 
starting from 100 seconds after the triggers with an on-board X-ray telescope (XRT) 
and a UV-optical telescope (UVOT). The mission will also issue prompt alerts to 
the ground to allow rapid follow-up observations from the ground-based telescopes. 

There are two important differences between the reverse shock emission and 
the forward shock emission. First, the reverse emission is prompt and short-lived. 
The electrons are continuously accelerated only until the reverse shock crosses the 
initial shell. This happens at the deceleration radius rdec, which is rr for the 
thin shell case, or at rx for the thick shell case (see H3.3(l . The observer's time 
at which this occurs is usually defined by the crossing time, tx — Tnax{tr,T), 
where tp = [{3E /AnTQnmpC^y^^ /2Tqc]{1 + z), and Fq is the initial Lorentz factor. 
The fireball Lorentz factor at the shock crossing time is Fx = min(ro,Fc), where 
F, ~ 125^5^ n-i/8T2~^/**[(l + z)/2]3/8. The reverse shock emission is therefore gen- 
erally divided into two segments. For t < tx, the synchrotron spectrum is (as in 
the forward shocks) characterized by a four-segment broken power law, as discussed 
in ij3.5l However, for t > tx, since no new electrons are accelerated, the emission 
above the cooling frequency {v > Vc) totally disappears. Second, before and at the 
crossing time, both the pressure and the internal energy density across the contact 
discontinuity are the same, but the particle density in the shocked shell is much 
larger than that in the shocked ISM. Given similar electron equipartition and injec- 
tion factors (ee and i^e), the typical energy per electron is much lower in the reverse 
shock region than in the forward shock region. So the synchrotron peak frequency 
in the reverse shock is much lower than that of the forward shock, mainly in the 
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Fig. 14. Typical GRB optical early afterglow lightcurve that includes the contributions from 
both the forward and the reverse shock emission components (from Ref.120). 



typically characterized by an optical flash and a radio flare. 

Compared with the forward shock emission, there are more physical cases to be 
considered in the reverse shock lightcurves, depending on issues such as whether 
one is in the thick shell regime (in which case the reverse shock becomes relativistic 
before i x ) or in the thin shell regime (in which case the reverse shock only becomes 
mildly relativistic at i x ); what is the relative position of the observational band with 
respect to the break frequencies both before and after the crossing time; fast or slow 
cooling; whether the medium is an ISM or wind-like, etc. For the ISM full 
discussion about various lightcurve cases (12 altogether) is presented in Ref. 11181 
A closer investiga tion taking reasonable parameters reveals that there are only 4 
most relevant case depending on whether the shell is thin or thick and another 
parameter TZ^, = vrJ Vm.ritx), which defines whether the observational band is 
above (for TZn > 1) or below (for TZi, < 1) the injection synchrotron frequency 
of the reverse shock emission at the shock crossing time. A common feature of 
all four cases is that the flnal temporal decay slope is steep, i.e. Fi, cx with 
a ^ 2. This segment of the lightcurve has been identified in GRB 99012o^ and 
GRB 02121 1'^^''^ , lending credence to the reverse shock scenario. The reverse shock 
lightcurve component eventually joins the forward shock component at a later time. 
A generic expected early optical afterglow lightcurve involves two peakJi^ (Fig. 
[T^ . i.e., a reverse shock peak with {tr,p, F,j^r,p) at the beginning of the a ~ 2 
lightcurve segment (which usually corresponds to tx), and a forward shock peak 
with (i/,p, F,^j,p), which is the transition point from the cx t^/^ lightcurve segment 
to cx t""'^ lightcurve segment corresponding to the typical forward shock synchrotron 
frequency crossing the band. 



optical/IR regimi 




and the reverse shock emission component is 
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The information provided by the reverse shoc k eniission has been widely u sed to 
estimate the initial Lorentz factor of the fireball Umzmimm 1 Usu- 
ally these methods rely on the poorly known shock parameters (cg, es and p) as 
constrained by the forward shock modeling, and they implicitly assume that same 
parameters also apply in the reverse shock. The absolute values of these parameters 
are usually used to derive Fq. We have recently proposed a straight forw ard recipe 
to derive Tq by using a combined reverse and forward shock analysi^^^. The only 
information needed is the time and flux level at both the reverse shock and the 
forward shock emission peaks. Since there ar e sirnple correlations of i^rm i^c and 
Fu,7n between both shocks at the crossing tim^^^^, the ratios of the two peak 
fluxes and the two peak times only depend on Fg (or a parameter depending on 
both Fq and Fc for the thick shell case), the unknown parameter TZu, as well as the 
ratios of the shock microphysics parameters. Suppose that both lightcurve peaks 
are detected in an idealized observational campaign; the two unknown parameters 
Fq and TZi, can be then solved for from the peak flux ratio and the peak time ratio, 
if the microphysics parameters are the same in both shocks (as implicitly assumed 
in almost all previous studies). Otherwise, the derived Fq (and 7?.^) can be still 
derived as a function of the ratios of the microphysics parameters. No knowledge 
about the absolute values of and es is needed. If the GRB central engine is 
strongly magnetized, it is natural to expect that the magnetic field in the reverse 
shock region could be stronger than that in the forward shock region. Generally, 
one can introduce a free parameter TZb = Br/Bf into the problem, which may be 
solved for along with Fq if enough information is available. Case studies suggest 
that the central engines of GRB 990123 and possibly of GRB 021211 are strongly 
magnetizec^^ (see also Ref. '143i) . This provides an independent clue in addition 
to the gamma-ray polarization^ (cf. Refs. |33l I44|l to suggest a magnetized central 
engine and fireball, although not necessarily Poynting-flux dominated. 

For a wind-l ike environment, the early afterglow lightcurves are considerably 
different! ' Due to the high density of the shell at the shock crossing time, 
the reverse shock emission is in the fast cooling regime. After shock crossing, the 
emission above Vc disappears. For typical parameters, Vc,r{t>i) is well below the 
optical band, so that the lightcurve after the crossing time is characterized by a 
steep decay (oc t~'^~^, (3 is the spectral index, which is either 1/2 or p /2 fo r the 
case we are interested in) due to the off-axis angular time delay effecll I This 
rapidly-decaying reverse shock emission joins the forward shock lightcurve later, 
at a time earlier than the one found in the ISM cas^^^H. There is also a simple 
correlation between the fluxes and times for the sh ock c rossing and the crossing of 
the forward-shock typical frequency across the bancPSI, which is relate d to F n and 
TZb 'AS well, so that a similar recipe to derive Fq and TZb as in the ISM cas can be 
utilized. When the self-absorption frequency is above injection frequency in the fast 
cooling regime (which is likely in some parameter regimes), balance between cooling 
and self-absorption heating implies a bu mp in the electron distribution spectrum, 
and hence, in the emission spectrunPa. This gives rise to interesting spectral and 
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temporal signatures for the early afterglows that can be used to diagnose fireball 
and wind par amet ers. Such a signature is also relevant for the case of a constant 
dense mediunA^. 

The significance of Fq and the forward to reverse field ratio TZb in understanding 
the GRB fireball content and prompt emission mechanism has been discussed in 
H4.ll and M.'M In the Swift era, an abundance of early afterglow data is expected, 
and it is desirable to systematically analyze these data to retrieve essential fireball 
parameters such as Tq and TZb- ^ follow-up statistical analysis of these data along 
with other measurable parameters would eventually lead to constra ints o n, or even 
the identification of, the GRB prompt emission site and mechanisnJ^^I 

The early afterglow data may also reveal dis tinct emission features from pro- 
cesses such as pair loading-3Sl and neutron deca}!^^. More detailed studies about 
these emission signatures (spectra, lightcurves) are n eeded in order to differenti- 
ate them from other signatures such as post-injectiorP^ and the passage of th e 
forward shock synchrotron injection frequency across the observational banci^i^. 



5.2. Short GRBs and other possible sub-categories 

In a discussion of prospects for future GRB studies, one major unsolved puzzle 
concerns the nature of the short, hard bursts, which comprise of about 1/4 of the 
total GRB population detected. Essentially all the information about (and from) 
afterglows discussed above is for the long bursts. The short bursts still remain 
as mysterious as the long bursts were before 1997. So far only upper limits on the 
afterglow emission of a few short bursts are available, alth ough the BATSE archival 
data contain marginal evidence of weak X-ray afterglow J '. Directly scaling the 
long burst afterglow models to the short bursts, the calculations are straightforward 
with most parameters unchanged except for a smaller total energy (due to the short 
duration if a roughly similar luminosity is assumed) , or a smaller energy per solid 
angle, if the jets are broader (as might be expected e.g. for NS-NS mergers where a 
large envelope is absent), and possibly a lower external medium density (as might 
be expected for the NS-NS merger scenario if these wander out of the galaxy) . The 
results indicate that the afterglows for short bursts are faint and consistent with 
the curren t upp er limits, and that the afterglows are most easily detected in the 
X-ray bancPSl 

On the theoretical side, simple estimates reveal that NS-NS mergers are likely 
to result in central engines with short durationJ^LSSO 153 Extensive numerical 
mo de ling: has been carried out to reveal how compact merger events can produce 

GRBjSMEnzEni. 

Significant coUimation is required for the neutrino-driven model 
in order to account for the detected isotropic-equivalent luminosity of short bursts 
(assuming cosmological distances) , while a magnetically-driven model requires fields 
as high as 10^^ G t o acc ount for the short duration (otherwise, mergers may also 
produce long bursts^P^. There is a suggestion that short bursts may also originate 
from collapsarJ^. Other ideas to account for the long and short bimodal distribu- 
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tion include whether or not strong GW losses are incurred in the spin down of a 
millisecond magnetar modCi^"^, and whether, in a coUapsar scenario, the black hole 
- torus involves hyper-accretion or suspended accretion because of the interplay 
between the disk and the black hole spin^Q-. All these models, however, have to be 
developed further to circumvent the apparent lack of the afterglows for the short 
bursts. 

Swift, thanks to its rapid slewing ability, is expected to be able to catch faint 
afterglows of some short bursts, if they exist, and consequently to lead to the 
identification of their locations, redshifts as well as other properties such as jet 
beaming. Merger events have been expected to occur in regions with a large offset 
from the host galaxy center (due to the asymmetric kicks during the formation 
of the NSs^^^; see, however, R.ef. I5()()|l . the position information alone relative to 
the host is already very essential for the short bursts' identity and/or merger rate. 
Furthermore, ambient density estimates through broadband afterglow modeling 
will also help to reveal whether short bursts are located in low-density regions as 
expected at least for some fraction of short bursts in the merger models. 

Will new sub-categories of GRBs be identified in the Swift era? One can only 
speculate on this. It may be unlikely to identify new categories based on one piece of 
information alone (e.g. burst duration). However, cross correlations among multi^; 
parameters may well lead to the identification of new c atego ries, such as XRFa^, 
long-lag bursts assoc iated with the super-galactic plan q^'" ! and the sub-energetic 
rapid-decaying burstJ^^. Alternatively, and equally exciting, it may be that new 
observational evidence and theoretical modeling might reveal that there is a unified 
paradigm behind all these apparently diverse sub-classes. 



5.3. High energy photon emission 

A new and so far barely explored window in the electromagnetic spectrum is the 
high energy extension in the GeV-TeV photon energy range, which holds signif- 
icant promise for a better understanding of GRB. After the better studied MeV 
prompt emission and the broad band low energy (X-ray down to radio) afterglow 
studies, the final (photon) frontier is the very high and ultra-high energy domain, 
extending at least up to TeV, and possibly beyond. In the intermediate range of 10- 
100 MeV observations, while not num erou s, have been mad e with relatively high 
significance, e.g. with SMMI^, OSSePHI, COMPTElPH, etc. In the GeV-TeV 
energy range, only modest to low significance detections of a handful of GRBs have 
been reported, e.g. with the EGRET spark chambeJ^] or TASC calorimeteiH2l or 
from ground-based air Cerenkov telescope^^^l^^. These results have wetted the 
appetite for obtaining higher quality data, and have led to a vigorous campaign 
of instrumental developments in h igh e nergy regimes. Several planned new mis- 
sions including A GILl P ^ , and several planned or already operating 
ground-based air and water Cer enkov telescopes (e.g. Milagro, VERITAS, HESS, 
MAGIC and CANGAROO-II, etcPSl), are expected to open a new era of unveiling 
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the mystery of GRB high energy emission. 

On the theoretical side, high energy photons in this energy range are expected 
from the leptonic component of the fireball (e.g. electron IC emission from various 
emission sites) as well as from the hadronic component of the fireball (e.g. proton 
synchrotron, 7r+ synchrotron emission and tt'^ decay from or pn,pp interactions 
in various emission sites (see H3.7.5I for a discussion in the external shock case). 
Within the standard fireball shock scenario sketched in Fig. 2, we could in principle 
have the following components emitting in the GeV range and possibly above. 



Electron self-IC component from the external forward shod* 
Electron self-IC components from the external rever se shqck or the cross-IC com- 
ponents between the reverse and the forward shock jza™. 

Electron self-IC component from the internal sho cks ^ '^[^'^^^ ^ 
Proton synchrotron emission in the external shock^Stl^^^^S 

Photo-meson cascade emission from the external 

shockJSSlHMl 

Proton synchrotron emission and photo-meson cascade emission from the internal 
shock&5Q9^; 

Cascade emission resulting from pn inelastic collisions during the early phase 
of fireball evolution when the neutron component decouples from the proton 
componenlP^^^^ 

Baryonic photosphere component (and possibly its Comptoniz at io'^coff P ouent ) 
extending into the GeV regime (for a low-cr high-entropy 



The above is not an exhaustive list, and not all of the above processes may be op- 
erative at any time or in all models. For example, if the GRB wind is strongly dom- 
inated by a Poynting flux so that the GRB prompt emission is not due to internal 
shocks, the internal-shock-related GeV components in the above list would be sup- 
pressed or absent. Within the coUapsar scenario, the precursor emission tends to be 
in low energies. But in the supranova scenario, the presence of the presumed super- 
nova shell and hot pulsar wind nebula provide copious photon or baryonic targets for 
additional IC and hadronic interaction components in the GeV rangJ^^ SBSSlffi Ol 
Such emission components may also exist for the type-II coUapsars such that the 
delay between the SN and GRB is short. 

An ideal theoretical approach would start with a study of all the above emission 
components within a unified theoretical GRB framework, to compare the relative 
importance of each component as well as the parameter regimes for each component 
to dominate. An example of the potential benefits of such an approach is provided 
by the recent observation of a distinct GeV component identified recently from an 
archival BATSEhurst, GRB 94101lE3, which evolves separately from the usual low- 
energy (sub-MeV) component. So far, only the high energy emission components 
in the external forward shock have been studied in a consistent mannei"^^^ 186] 
In the external shock scenario, the proton synchrotron and hadronic cascade emis- 
sion components are only important when es is close to equipartition, while eg is 
very small, i.e. or order of me/mp or below. This is especially so when the proton 
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Fig. 15. Model prediction of GRB high energy (GeV) afterglows as compared with the GLAST 
sensitivity. For the IC-dominated parameter regime, an extended (hour-long) GeV afterglow should 
be detectable by GLAST for bursts at typical cosmological distances (from Ref.186). 



power-law index is taken a s 2.2 -2.3 rather than the special value 2 adopted in most 
of previous 

studie^22L322|MSI. 

For the more commonly invoked parameters (e.g. 
Ee ^ 0.1, eB ~ 0.01) as inferred from afterglow studies -'^'^-'^ l'S2^ hadronic com- 
ponents are completely buried beneath the electron IC component, which forms a 
distinct bump in the GeV regime. The MeV-GeV lightcurve predicted wi thin this 
parameter regime reveals a second broad bump lasting hours or even a d a\ l^^^ l ^ ^, 
and the flux level is detectable by GLAST for typical bursts at z ~ pjM (Fig. 
I15|l . According to this scenario, the previously detected long-lived GeV emission 
for GRB 9402 llEl 

is simply a nearby burst whose GeV afterglow is bright enough 
to be ca ught by EG R E 2 '^"^ ^8(>^ although alternative interpretations have been also 
proposecl"^^^^^^^. GLAST would reveal whether such kind of long-term GeV 
emission is common, which would pose important constraints on GRB shock phys- 
ical parameters °". 

Another prominent issue is the high energy photon cut-off in the GRB spectra. 
During the escape of high energy photons, they are subject to absorption through 
77 interactions with low energy photons within the source and in space. These 
interactions produce electron-positron pairs and greatly degrade the photon flux 
level in the original energy range received at the detector. At the source, such 77 
pair process would result in a well-defined high energy cutofl " that would serve 
as an important diagnostic of the fireball initial Lorentz factoi'^S-'^. Another con- 
sequence of this process is that the produced pairs are expected to synchrotron 
radiate again within the local magnetic fields. Depending on the compactness of 
the emission region, the secondary emission process may modulate the emergent 
spectrum u nder cer tain conditions and possibly smear out the initial distinct spec- 
tral featureJ^SI^. A self-consistent pair-modulated numerical model for the GRB 
prompt emission is needed for this. Even if a high energy photon is not absorbed, 
it may be scattered by the resulting pairs and even by the electrons associated 
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with the fireball baryonJ^MI Furthermore, if the source is not compact enough so 
that the T eV p hotons can escape, they may be still absorbed by the cosmic IR 
backgroun d^ ' ' '^ 1 The mean free path of the TeV photon depends on both the GRB 
and the IR background models. Typically a TeV source cannot be detected beyond 
z ^ 0.1^1^ so that the Milagrito event for GRB 970417-"^^^, if real, has to come 
from a nearby source whose compactness is not high (e.g. from external shocks, or 
with a high Lorentz factor). For typical cosmological sources, TeV emission will be 
almost completely absorbed by the IR background. The resultant pairs, if within 
a not-too-strong intergalactic medium (IGM) magnetic field, would IC up-scatter 
the cosmi c microwave background photons to produce a delayed, GeV emission 
componenJSnESEni. A well-measured spectrum may be used to infer the poorly 
known strength of the IGM magnetic field (see Ref. i516i for a detailed modeling 
within the context of blazars). 

5.4. Ultrahigh energy cosmic rays 

Two important non-electromagnetic channels in which GRB may be prominent 
sources are cosmic rays and neutrinos, neither of which have so far been measured. 
GRB models involving shocks as sites to accelerate electrons which produce prompt 
gamma-rays and long-term afterglows naturally suggests that baryons (most likely 
protons) should be accelerated by the same shocks as well. These accelerated ions, 
if not bound in the system and not destroyed during their propagation, would arrive 
at Earth and be detected as cosmic rays. Those trapped within the system would 
interact with photons and other baryons to produce high energy neutrinos that 
might also be detected from Earth. 

Discussions about GRBs as cosmic ray accelerators are mainly aimed at explain- 
ing the portion of the cosmic ray spectrum above 10^^ eV, the "ankle" region where 
the spectrum starts to harden, and in particular the ultrahigh energy cosmic rays 
(UHECRs) above 10^" eV, also referred to as Greisen-Zatsepin-Kuzmin ("GZK") 
cosmic rays. The isotropic distribution of their arrival directions and small magnetic 
deflection they suffer at these high energies suggest their extra-galactic origin, and 
the requirement that they must survive the attenuation by the cosmic microwave 
background through photomeson interaction const rains the ir distances to radius 
of about 50-100 Mpc, the so-called "GZK" volumJSnEIHI. 

Two broad classes of 

models have been suggested to interpret the UHECRs, the "top-down" scenarios 
that attribute UHECRs to decay products of fossil Grand Unification defects, and 
the "bottom-up" scenarios that suggest UHECRs are hadrons accelerated in astro- 
physical objects to these high energies. Among the few serious viable candidates 
in the bottom-up scenario, GRBs (and/or perhaps AGNs) arc considered to be 
likely sources which may fulfill the constraints on the known data-*^^^ 169 ol9^ Two 
sub-scenarios of the bottom-up GRB cos mic ray origin have been suggested, one 
involving acceleration in the internal shockj ^^l'^^l'^ and one involving accelera- 
tion in the external shoclllS2»^22j The suggestion is based on two coincidences, i.e., 
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the shock conditions required to accelerate protons to ^ 10^° eV are similar to the 
conditions required for generating the observed prompt gamma-rays, and the ob- 
served UHECR energy injection rate into the uni verse (--^ 3 x 10'*'* erg Mpc^'^ yr^*) 
is similar to the local GRB energy injection rat^^^^^. Both coincidences have 
been questioned on various grounds ''^'^ ^'^'^ 5 JSE2EI^ but these have been met by 
effective counter- argu ments, using new data and additional considerations for both 
the internal shoclJSniand the external shock"^ scenarios. At this time, GRBs are, 
and remain, a promising candidate for the UHECR origin. However, both the inter- 
nal shock and the external shock scenarios have major caveats, which have nothing 
to do with their ability to accelerate cosmic rays but rather with the generic shock 
itself. The internal shock scenario relies on the assumption that the GRB prompt 
emission is due to synchrotron from electrons accelerated in internal shocks. Al- 
though this is the leading scenario, there is no direct proof so far, unlike in the case 
of the external shock, where this origin of the radiation is quite convincing. For 
example, a Poynting-flux-dominated GRB model would be incompatible with an 
internal shock origin of both synchrotron 7-rays and UHECRs. The extern al sho ck 
model, on the other hand, may h ave t o rely on a magnetized ex tern al mediunJS22 (as 
expected in pulsar wind bubble^ED in the supranova scenaric^l) in order to reach 
the desired cosmic ray energy. On the other hand, the radiation and acceleration 
responsible for 7-rays is itself unexplored in Foynting-dominated scenarios, and the 
need for such a scenario is at the moment not proven, while the supranova scenario 
is incompatible with the observations of GRB030329/SN2003dh. The uncertainties 
surrounding possible alternative scenarios are significant, and the standard shock 
scenario remains the most amenable to quantitative modeling and testing, which 
should be able to provide useful constraints on its ability to explain the data. 

A direct proof of the GRB-origin of UHECRs is not easy. The next generation of 
cosmic ray detectors, such as the Auger Observatory, will have a substantially en- 
hanced effective target area, which will greatly improve the cosmic ray count statis- 
tics. It also combines elements of the two techniques which currently lead to different 
results, namely air fluorescence telescopes and water Cerenkov surface tanks. This 
will help to disentangle the two distinct scenarios (top-down or bottom-up) and 
to reveal whether a GZK feature indeed exists. Within the bottom-up scenario, 
the direction information may prove or significantly constrain the AGN model (the 
close competitor of the GRB model) and eventually shed light on whether GRBs 
are indeed the sources of UHECRs. 

5.5. High energy neutrinos 

Regardless of whether GRBs can accelerate protons to ultrahigh (GZK) energies, 
they must be able to accelerate protons to some high energies. The implication is 
that high energy neutrinos and high energy photons (as discussed in ^5.3|l must ac- 
company the current emission seen at sub-MeV energies. Widely discussed processes 
for high energy neutrino emission include 
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• p7 process: P7 — > A+ n7r+ — > ne^i'f.i^^i'^] 

• pp process: pp tt^/K"^ . . . ^ fiiy/^, . . . — > eVeV,_i.Vfj. . . .; 

• pn process: pn —f / /xf^ . . . — > ez^gp^i/^ . . . 

The dominant p7 process occurs at the A-resonance, which has the threshold con- 
dition epe^>,0.3 GeV^ in the center of mass frame. In the case of GRBs, this is 
usually translated in the observer's frame to £^6-^ >,0.3 GeV^F^ when both protons 
and photons are generated in relativistic shocks (e.g. in the internal or external 
shock scenarios), or to epe^>,0.3 GeV^F when the photons are generated in SN 
shells while protons are accelerated from relativistic shocks (e.g. in the supranova 
model), where F is the bulk Lorentz factor. The threshold condition for both pp 
and pn interactions is that the relative drift energy between these baryons exceed 
the pion rest mass, i.e., e' > 140 MeV. Since both p and n have a rest mass close 
to 1 GeV, the threshold of pp and pn interaction only demands semi-relativistic 
relative motions. 

In a GRB event, there are multiple sites where neutrinos with different ener- 
gies are generated. Below is an non-exhaustive list which encompasses most of the 
processes discussed in the literature, in a sequence of ascending neutrino energy: 

• MeV neutrinos: If GRBs are originated from stellar collapses, they should be 
associated with strong thermal MeV neutrinos like supernovae. In most models, 
the collapse results in a black hole - torus system, and the thermal neutrino anni- 
hilation is one of the leading processes for launching the fireball. However, these 
thermal neutrinos are extremely difficult to detect from cosmological distances, 
due to the very low cross section for i/N interactions at these energies (the cross 
section increases cx at lower energies, and oc E^, at higher en ergies) . 

• multi-GeV neutrinos: GRB fireballs may be neutron-riclPSESSI During the fire- 
ball acceleration phase, neutrons can decouple from protons when the elastic 
scattering condition breaks down. The relative drift between b oth s pecies re- 
sults in inelastic pn interactions giving rise to 5-10 GeV neu trino J^^^ I A similar 
process also occurs within sub-photospheric internal shockJ^^, which extends 
significantly the parameter space for the inelastic neutrino collision condition. It 
was suggested earlier that pp interactions within the internal shocks can also 
give rise to a 30 GeV neutrino burst, although magnetic fields can inhibit the 
inter-penetration of charged species streams with different velocities; 

• multi-TeV neutrinos: Within the coUapsar scenario, the relativistic jet launched 
from the base of the flow (presumably the black hole and the torus) has to 
penetrate through the stellar envelope before breaking out and generating the 
GRB. The internal shocks below the envelope accelerate protons that interact 
with thermal photons within the envelope (i.e. pj interaction). Regardless of 
whether the jet finally penetrates through the envelope or gets choked, it will 
generate strong multi-TeV neutrino sig nalM The signature is enhanced or 
even dominated by pn,pp interactions and could be used as a diagnostic about 
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Fig. 16. The diffuse muon neutrino flux for both the coUapsar and the supranova models (from 
Refs.463,527). 



the type of progenitor starJS^ 

• PeV neutrinos: P7 interactions within the conventional internal shoc ks which 
produce prompt gamma-rays typically generate 10^** - 10^^ eV neutrinoJEOESHl 

• EeV neutrinos: P7 interactions within the external reverse shock give rise to even 
higher energy neu trinos. For a constant density medium the typical energy is 
10^^ - 10^^ e\E2ni, while for a wind med ium , the typical energy is 3 x (10 — 
10^^) eV and extending above i j530|53j | jj^ ^j^g forward shock region, assuming 
the blast wave can accelerate protons to ultrahigh energie j ^ ^^'^'^^^ a neutrino 
afterglow is expected with the peak energy ~ 10^^ g\/:531 532 

In the supranova progenitor scenario, the existence of the pre-supernova shell 
provides extra targets for additional neutrino signals. These signals are broad-band 
with a distinct 1-10 TeV fea ture with high co unt rates, which is a useful clue to 
test the supranova 

hypothesiJSSnniEnaiMi. 

The neutrino model predictions for 
GRB 030329 indicate that a 2-day delay supranova model predicts a one order of 
magnitude higher neutrino event rate in the TeV-PeV regime than the conventional 
burst, a signal which is detectable by the full ICECUBE as an individual source 
at the observed redshift z = 0.17. If it were to be detected by under-ice neutrino 
detectors at the South Pole or underwater detectors such as An t ares , this would 
further significantly constrain the validity of the supranova model^^. 

All the above are theoretical predictions. So far there is no report about the 
solid detection of neutrino signals from GRBs. Current and future ice or water 
Cerenkov neutrino detectors, AMANDA-II, ANTARES, and ICECUBE will make 
it possible to detect GRBs from this new channel. If detected, neutrinos from GRBs 
could be also used to study neutrino physics itself as well as to have ramifications 
for cosmologjl^i. 



5.6. Gravitational waves 

Binary compact-object mergers, such as NS-NS, NS-BH, BH-BH mergers, as well as 

BH-WD, BH-Helium star etc.. have lo ng been con sidered as possible sources o f gray- 

itational waves (GW )535 | 536 | 537 | 173 | 350.538 | T71| ^ ^ave core-collapse event J&^'^ l ''^40 | 541 | 174 | 542 | 543 | 544 | 
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Fig. 17. The gravitational wave strains for tiie merging double neutron star model and the 
coUapsar model as compared with the sensitivity of the advanced LIGO (from Ref. 174). 



This has been regardless of whether they could produce GRBs. However, since these 
events are also the leading candidates for being GRB progenitors, it is natural to 
expect that a GRB is associated with a GW burst (although the contrary may 
not hold). For a selection of GRB-related GW literatures, see e.g. Ref. 15451 A co- 
incidence between a GW signal and a gamma-ray signal would greatly enhance 
the statistical significance of the formei^^^, making searching for GRB-associated 
GW bursts of great interest. A binary coalesce nce proc ess can be divided into 
three phases, i.e., in-spiral, merger, and ring- For coUapsars, a rapidly 

rotating core could lead to development of a bar and to fragmentation instabil- 
ities that would resemble similar GW signals as in the binary merger scenarios, 
although a larger uncertainty is involved since the fragment masses and coherence 
times are unknown. The GW frequencies of the various phases cover the 10— 10'^ Hz 
band which is relevant for the Laser Interferometer Gravitational-wave Observ atory 
(LIGO) and other related detectors such as VIRGO, GEO600 and TAMASO^^. 
Because of the faint nature of the typical GW strain, only nearby sources (e.g. within 
~ 200 Mpc for NS-NS and NS-BH mergers, and within ~ 30 Mpc for coUapsars jEl 
have strong enough s ignals to be detectable by LIGO. When event rates are taken 
into accounlllSniMni^ 

order of magnitude estimates indicate that after one-year of 
operation of the advanced LIGO, one in-spiral chirp event from a NS-NS or NS- 
BH mer ger, and probably one coUapsar event (subject to uncertainties) would be 
detectecffZaESI. Other binary merger scena rios such as BH-WD and BH-Helium 
star mergers are unlikely to be detectabl^^, and they are also unfavored as sources 
of GRBs according to other argument&P*^. 

The detection of a GRB-GW burst association would have profound implications 
on GRB studies, and could shed light on several unsettled issues as discussed in 21 
First, it would help to identify the GRB progenitor. GW signals in the advanced 
LIGO band are different for the merger and the core collapse scenarios, so that a 
coincident GRB-GW detection would unambiguously differentiate both scenarios. 
Furthermore, two GW bursts are expected for the supranova model, with the second 
one coincident with the GRB. Either detection or non-detection of the precedent 
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GW signal would provide additional information on the validity of the supranova 
model (in addition to the negative evidence provided by GRB030329/SN2003dh). 
Second, it could help to settle the debate about the location of GRB prompt emis- 
sion (i.e. external or internal, M.2|l . Because the GW wave is associated with the for- 
mation of the central engine (and therefore the launch of the GRB jet), a short time 
delay (<,0.1 s) for the GRB emission with respect to the GW emission would favor 
an internal scenario (e.g. the internal shock model), while a longer delay 10— 100 
s) would favor an external model (e.g. the external shock model) ^^'^ (this how- 
ever ignores the delay incurred by the jet in punching though the stellar envelope, 
which could itself be of order 10-lOOs); Third, it would help to pin down the GRB 
jet configuration. It has been suggested that the GW polarization information could 
be used to infer geometric information such as the angle betwe en th e rotation axis 
(presumably the center of the jet) and the line-of-sig ht directioJ^^. Such informa- 
tion, when combined with other information from prompt gamma-ray and afterglow, 
could help to confirm whether and how the GRB jets are structurcQ'*^^' ^ 136 . 137 | 



5.7. GRBs & cosmology 

One of the most exciting prospects for future GRB studies is that they could become 
a unique tool to investigate the high redshift universe. Although this possibility is 
currently mainly a theoretical expectation, the prospects look bright. The optimism 
is underpinned by two aspects of the evidence obtained so far. 

Clearly, since GRBs are stellar events, some GRBs are expected to exist at 
(much) high er re dshifts than z = 6.4, the current high redshift record held by a 
SDSS quasaJS^. This is supported by several lines of reasoning. (1) Preliminary 
polarization data on the cosmic microwave background collected by Wilkinson Mi- 
crowave Background Probe (WMAP) indicate a high electron scattering optical 
depth, hinting that the first stellar objects in the universe should have formed 
as early as z ~ 20'''''^; (2) Independent theoretical simulations of the formation 
of the first stars sim ilarlv conclude that these should have formed at redshifts 
(15-40 jSS4|555l. 

Because there is convincing evidence that at least long GRBs 
are associated with the deaths of the massive stars, it is conceivable that high-z 
GRBs {z<lh — 20 or even higher) exist. In fact, although subject to considerable 
model uncertainties, theoretical modeling of GRB redshift distribution based on 
the standard ACDM universe suggests that >,50% of all GRBs on the sky originate 
at z>,^^3, (3) Using empirical luminosity laws such as gamma-ray variability or 
spectral lag correlations with fiuence, rough redshifts can be derived for a large 
sample of the BATSE bursts, and it is found that a good frac tion of these bursts 
have redshifts in excess of z ~ 6, and some even in excess 

The emission properties of GRBs and afterglows offer several unique advantages 
for studying the high-z universe (see als o Ref. [ 5 5611. (1) The luminosities of the 
prompt gamma-ra^l ^ and the afterglow^^^^fdo not fade rapidly with increas- 
ing redshifts (unlike for quasars). In fact, the highest-redshift burst identified so far 
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(GRB 000131 at z — 4.5) does not show a substantially fainter luminosity in both 
the prompt emission and the afterglow. This is due to a favorable combination of the 
cosmological time-dilation and the flux decay resulting in a positive K-correction 
effect. For afterglows, the high-2 bursts are especially favorable for d etect ion in the 
IR band, thanks to the bright reverse shock emission at early epoch^^El (2) GRB 
afterglows have intrinsically high luminosities (10^^ — 10^'* erg s~^). Even at moder- 
ate redshifts, they already greatly out-shine the host galaxies. This would especially 
be the case at higher redshifts, as the host galaxies are expected to become even 
fainter; (3) The intrinsic power-law spectrum of GRB afterglows greatly simplifies 
the extraction of the IGM absorption features. 

The implications of high-z GRB studies for cosmology include the following. 
(1) WMAP only provides an integral constraint on the reionization history of the 
universe. Since both the number density and the intrinsic luminosity of quasars are 
expected to fade rapidly beyond z ~ 6, only GRBs and afterglows may be able 
to act as bright beacons to illuminate the end of cosm ic da rk age , and 

to probe the reionization history of the early universJ^^ (2) Although the ap- 
parent GRB luminosity is by no means a standard candle, the beaming-corrected 
GRB luminosity appears to be standar j^''-*^^-*^'^^^'^"^^* . Also empirical luminosity 
laws^*^^ ^-'^^j if confirmed, could be used as Cepheid-like correlations. These raise 
the possibility of using GRBs to derive a Hubble-diagram and to constrain the cos- 
mological parameterj2SSIlM01 although the current data are too sparse to draw a 
firm conclusio n: (3) Pr eliminarv evidence on the evolution of GRB properties has 
been sugg If this is confirmed by future photometric redshift measure- 

ments of a large amount of GRBs (with selection effects properly taken care of), 
the data could be used to infer the possible redshift-evolution of the GRB pro- 
genitors; (4) Future extensive afterglow monitoring for many bursts would help to 
constrain the local environments of GRBs as well as their redshift evolution. In 
particular, it would tell whether GRB afterglows are decelerated by the IGM (with 
an increasingly high density at higher reds hifts) or by a stratified constant medium 
bubble cleared by the progenitor stai^^ (5) Insights into cosmological structure 
formation and star formi ng hi storv would be gained through studying distributions 
of the GRB host galaxieJSSl. 

A great challenge in the Swift era is to get direct redshift information for high- 
z bursts. The most straightforward approach is to search for the Lya break (at 
1.216/im[(l + z)/10]) and the Gunn-Peterson trougll^^in the IR band. However, 
since the Swift UVOT wave-band extends only to ~ 0.6/im in the long wavelength 
regime, the redshifts of GRBs with z > 5 can not be directly identified with Swift. 
This poses however a great opportunity for ground-based IR cameras to follow up 
Swift-triggered GRBs within a short period of time. Several IR follow-up teams have 
developed such instruments and plan to carry out Swift follow-up observations in 
a timely manner. Besides IR follow-up, it is also pos sible that X-ray lines could 
become important distance indicators at high redshift J^^. Alternatively, redshifts 
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for GRBs whose ra dio after glows are unusually bright may be identified with 21 cm 
absorption featureJ^^^E^. 

Finally, using cosmological GRB high energy data, one may constrain a number 
of ideas in fundamental phys ics, e.g. related to quantum gravity effects, Lorentz 
invariance violations, etc! ^ ' '^'^ ! 

6. Summary 

Our understanding of GRBs has been greatly advanced since their discovery about 
30 years ago. Extensive observational efforts have revealed a rich phenomenology of 
both prompt emission and afterglows O, and a successful theoretical framework 
has been set up, which is able to interpret most of the observational data so far 
(SjSl- However, there remain many questions at the present stage of GRB studies, 
which greatly stimulate further observational and theoretical efforts O. Some 
of these questions will undoubtedly be addressed in the coming years with the 
advent of a number of upcoming space- and ground-based experiments, both in 
the electromagnetic and non-electromagnetic channels (^. As suggested by past 
experience, new challenges and surprises are bound to emerge, which will stimulate 
further extensive observational campaigns and theoretical efforts. The GRB field is 
likely to remain one of the most active fields in contemporary astrophysics in the 
next decades. 
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AST0098416 and the MoncU Foundation. We thank S. Kobayashi, T. Lu, D. A. 
Frail, J. P. Norris, E. D. Feigelson, Z. G. Dai, Y. F. Huang, and S. Razzaque for 
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